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Magnetic Susceptibility of Ferromagnetic Minerals 
Contained in Igneous Rocks. * 


By Syun-iti AKIMOTO 
Geophysical Institute, Tokyo University 


Abstract 


The thermal change in susceptibility of ferromagnetic mineral 
grains separated from igneous rocks was measured by a ballistic 
method in a weak magnetic field. The thermal change of susceptibility 
was classified into four types. These types were discussed in connec- 
tion with the condition under which the original rock was formed. 
The dependency of the Curie-point and the magnitude of the suscep- 
tibility upon the chemical constitution was examined, the content of 
TiO, in ferromagnetic mineral being especially taken into account. 
Besides, the relation between the magnitude of susceptibility and the 
grain size was obtained from the viewpoint of the magnetic property 
of small particles. 
§ 1. Introduction 
The magnetic properties of igneons rocks have been examined fairly in detail 
by a number of investigators (1), (2), (3), (4), (5), (6), especially in connection with the 
local geomagnetic anomalies and their changes. The magnetic properties of rocks, 
however, are due to those of ferromagnetic minerals contained in them, the other 
minerals scarcely contributing to the magnetism of rocks. In other words, the magne- 
tism of rocks is composed of the ensemble of magnetization of smail grains of rock- 
forming ferromagnetic minerals. So that the fundamental problem in rock-magnetism 
will be the general description of rock-forming ferromagnetic minerals from the 
physical standpoint in reference to their chemical and mineralogical characteristics. 
In the present study, the change with the temperature in magnetic susceptibility 
of rock-forming ferromagnetic minerals, which were separated from various kinds of 
mother igneous rocks, was measured, and the results were analysed in relation to the 
chemical composition, grain size and mineralogical characteristics of these grains of 
minerals. The aims of the present study are chiefly the following two. 
(1) The natural ferromagnetic minerals contained in igneous rocks are generally 


of chemicel constitution different from a pure stoicheometric magnetite. It will be 
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significant, therefore, to find the relations of the Curie: point and the magnitude of 


susceptibility to the chemical composition. 
(2) On the other hand, it has been expected that the dimensions of grains play 


‘an important role upon the magnetic properties such as initial permeability and 


* Contribution from Division of Geomagnetism and Geoelectricity, Cbesiveical Institute, Tokyo 
University. Series I]. No. 20. 
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Table 1 


No Tr toler nai Petrological Descripion Grain Size} Fe,0;+FeO TiO2 
1 | Niizima Phenocryst magnetite in biotite 
liparie 410p. | 63.01 31.46% |4.06 (*4.06)% 

2 |Lava in Asi- | Phenocryst magnetite in hyp. 

gara layer horn. dacite ' batou) 67.60 10.81 (13.73) 
3 Social vble: Groundmass magnesite in oli- 

Ito vine basalt 53 56.66 14.44 (20.08) 
4 | Simohutago- Phenocryst magnetite in hyp. ; 

yama, Hakone | aug. andesite 120 75.80 10.86 (12.46) 
5 | Kaziya Yuga- | Phenocryst magnetite in hyp. 

wara ' | dacite hottest Aes 74.44 10.24 (12.08) 


6 | Yugawara Phenocryst magnetite in aug.. 3 : 
hyp. andesite 130 69.84 8.73 (11.05) 
7 | Karataki, Phenocryst magnetite in aug. 
Hakone _| hyp. andesite 110 of 12.17 (13.79) 
8 | Hiroawara Phenocryst magnetite in aug. 
bearing hyp. dacite 130 83.28 11.40 (12.02) 
9 | Odawara Phenocryst magnetite in a 
| pumice 280 88.13 6.43( 6.79) 
10 | Taga vole. Phenocryst magnetite in oliv. y : 
aug. hyp. andesite 200°- | 75.60 10.73 (12.34) 
11 | Haruna volc. | Phenocryst magnetite in 
hyp. horn. dacite pumice 390 53.15 BL.93-c|. Gseatereass 
12. | Sukumogawa_ | Groundmass magnetite in and- ; ; 
Hakone esite 43 65.05 16.59 (20.18) 
13 | Okata Oosima | Groundmass ese in ; bicep tse 
a tholeiite | ; 63 47.19 33.71 17.20 (17.20) 
14 | Taga vole. Groundmass ssegneti in : : . 8 
olivine basalt 46 36.03 43.89 | 18.02(18.02) — 


Taga vole. Groundmass magnetite in 


oliv. aug. hyp. andesite 29, Sieis = A901 | 9.26 (15.79) 


tT * “This value 1 was calculated as 100% in total. — 
at ioe HK ed, 
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this table, the ferromagnetic minerals concern- 


Out @us Ker—w°0>| Xo ed here are the ferromagnetic ferrites chiefly 
ol composed of FeO, Fe.O; and TiOs. 

575°C 9.66 9.66 The petrological research of the original 

565 7.41 8.21 igneous rocks has been carried out by H. Kuno 


chiefly from a standpoint of analysis of pyro- 
455 Se0G te ae 72 3.28 xenes contained in them. (9) According to his 
opinion, No. 1 and No. 2 specimens in Table 1 are 


570 | 520 4.27 | 4.47 , 
the magnetite-grains in those rocks which were 
600 390 5.18 5.70 cooled slowly from a low temperature, being 
530 440 6.60 7.31 kept in the complete equilibrium, while No. 12, 
No. 13, No. 14 and No. 15 specimens are the 
eg =e ges magnetite-grains in those rocks which were 
540 450 7.84 8.70 cooled rapidly from a relatively high tempera- 
ture and No. 4, No. 5, No. 6, No. 7, No. 8, No, 9 
580 | 450 4.22 | 4.35 a 
No. 10 and No. 11 specimens are the magnetite- 
550 460 | 5.44 5.80 grains in those rocks, the.character of which 
570 7.72 _ lies between the above two cases. 


§ 3. Apparatus 


Be0 a6 The apparatus for measuring the thermal 

+580 1.54 change of the magnetic susceptibility of ferro- 
oe SET magnetic mineral grains in a weak field is the — 

— 600 1.67 same in its principle as that used: in ‘the pre- 

'. 550 ' fy 1.74 : | 9 24 vious studies upon : igueous rock samples. (3); (4), 


: (6) ie. the maguetization of specimen was meas- 

aa By. a ballistic ‘fiethod: at any temperature ‘in a weak magnetic field a few times 
' as much as that of geomagnetic field. 

ete est specimens were packed into a silica tube 3mm in inner diameter and 

. in : In the case of measuring the magnetization of the te aaa 

ils “the demagneting field is practically negligible in “comparison with the 

id. The electric ttc is -a ute "and non-magnetic 
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Fig. 1 No. 1 H=1.13 Oe Fig. 1 No. 2 H=1.35 Oe 


in this figure, the mode of thermal change of susceptibility of ferromagnetic minerals _ 


is similar so much to the previous results in the cases of the igneous rocks. It seems 
that there are four types of the mode of the thermal change of susceptibility. The first 
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Fig. 1 No. 14 H=1.35 Oe , Fig. 1 No. 15 H=1.58 Oe 


process it increases gradually up to about 450°C and then decreases abruptly until 
the ferromagnetism disappears, and the change in cooling process follows almost the 
same curve reversibly. This type of change will be called here the “ high-suscep- 
tibility and reversible type.’ The specimens belonging to this type are No. 1 
and No. 2 specimens in Fig. 1. The second is such the type as given by the X—T 
curve of the specimens No. 3, No. 4, No. 5, No. 6, No. 7, No. 8, No. 9, No. 10 and No. 11 
in Fig- 1, in which the magnitude of susceptibility has the same order of its value 
as the first type and the change in heating process is similar to that of the first type 
but the cooling curve differs distinctly from the heating one, and besides in re-heat- 
ing and re-cooling processes the susceptibility changes reversibly along the same 
curve as that in the previous cooling curve. This type of change will be named the 
“high-susceptibility and irreversible type.” The third is the type’ in which the 
magnitude of susceptibility is very small in comparison with the above-mentioned 
two types and it changes with temperature reversibly. This type of change will be 
named the “low-susceptibility and reversible type.” The specimens No. 12 and No. 13 
in Fig. 1 seem to belong to this type. The forth type is that the magnitude of 
susceptibility has the same order of its value as the third type and it changes with 
temperature irreversibly. This type of change will be named the “ low- susceptibility 
and irreversible type.” The specimens No. 14 and No. 15 in Fig. 1 seem to belong 
to this type. Thus, it will be summarized that there are two types “high-suscep- 
tibility ” and “low-susceptibility ” with respect to the magnitude of susceptibility and 
each of them is classified into two types “ reversible” and “irreversible” with Tes- 
pect to the thermal change of susceptibility. 

Comparing the above-mentioned classification of the %— T curves of the 
specimens with their petrological classification made by H. Kuno, it will be noticed 
that the specimens of “ high-susceptibility and reversible type” and of “ low-suscep- 
tibility type” of X—T curves correspond exactly to those cooled slowly from a low 
temperature and to those cooled rapidly from a high temperature respectively at the 
stage of the formation of their mother rocks, while the specimens of “high- -suscep: 
tibility and irreversible type” "correspond to those of petrologically intermediate tyne 


mentioned in &- 2. 
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It seems that the correspondence mentioned above can be physically explained 
by taking into consideration the characteristics of the Curie-point and the magnitude 
of susceptibility of these mineral grains dealt with in the forthcoming paragraphs. 
§ 5. The Curie-points 

The Curie-point, 9., of the specimens was defined to be the temperature at 
which ferromagnetism disappears apparently in the —T curve. As for the specimens 
of the irreversible type, there are two Curie-points, 9,1, and 0.5, corresponding to 
the heating and cooling curves respectively, where 9. is always higher than 6.2 in 
the present case. The values of the Curie-points thus determined are tabulated in 
Table 1. 

The Curie-point is a physical constant of each ferromagnetic material. The 
Curie-point of pure magnetite, for example, has been determined to be at about 580°C. 
When FeO or TiO, makes a solid solution with FesO,, some amounts of decrease in 
the Curie-point will be observed. In order to examine the relation between the 

: Curie-point and the chemical constitu- 
a tion of the specimens, 9, of fifteen 
600 specimens were plotted against the con- 

(0) tent of TiO. in them in Fig. 2. In the 
$00 Me case of the irreversible type specimens, 
0.2 was taken as 9,, since it may be 
on O the more stable and consequently sig- 
nificant Curie-point. It will be seen in 
Fig. 2 that the Curie-point of the four 


300 ae Pipes : , 
5 70 Ti 20 Ti Oh reversible type specimens is always 


Fig. 2 around 580°C and does not change with 
© reversible type specimen the content of TiO., while the 9.2 of 
O_ irreversible type specimen the irreversible type specimens de- 

creases as the content of TiO, increases. - | | 
The typical each two specimens among the two groups of the specimens of 
“reversible” and “irreversible”. types were chemically separated into ferric and 
ferrous oxides and titan oxide. They are No. 1 and No. 13 in the “ reversible” type 
specimens, and No. 11 and No. 14 in the “irreversible” type ones, the contents of 
FeO, Fe,0; and TiO, in them being given in Table 1. In the two typical specimens 
of the “reversible” type, these three kinds of metalic oxides can be interpreted to 
compose FeO. Fe.0; (magnetite) and FeO. TiO, (ilmenite) with very little excess. On 
the other hand, those of the two specimens of the “irreversible ” type Show a_ fair 


amount of FeO as an excess in ee above- mentioned calculation of Norm of magnetite 
and ilmenite. 


It may be assumed that the ferromagnetism of the 5 peversiule ” type specimens 


is practically. that of magnetite alone, because the magnetism of. ilmenite, which may 
simply be mixed with magnetite in this case, is far feebler compared with that of 
magnetite, tht Curie- -point being kept around 580°C regardless of the content of ilmenite. 


Fig. 3 shows X—T curve of an artificial V 
magnetite, where the Curie-point | is : 
580°C, and the X—T curve is rever- 
sible with respect to change in tem- 
perature. This *—T curve can be 
considered to be the standard curve of at! 
small grain sample of almost pure 
magnetite. In the “irreversible” type 
specimens, on the other hand, FeO, Ae 
FeO; and TiO: may compose complex 
ferromagnetic minerals considerably noe 
different from a pure magnetite and 
consequently the Curie-point decreases 
with the increase in the apparent con- 
tent of TiOs. ed 
§ 6. Magnitude of susceptility 

It will be noticed that the speci- 
mens given in Fig. 1 have a wide 
variety of the magnitude of % at room 
temperature. For the purpose of 


0 200 400 (obo 
Fig. 3 H=1.35 Oe 


studying the effects of impurity and grain size upon the magnitude of %, the values 
of X at the room temperature were examined in connection with the content of TiO: 


Grain sige and grain size, the result being shown 
=a in Fig. 4. That the magnitude of x 

decreases as the grains become finer and 

ey ag er as the content of TiO, increases will be 

fairly conceivable from the result given 

306 ee in this figure. However, there is also a 

close relation between the grain size 

200 ao and the TiO: content as will be seen in 

oan! *: Fig. 4. Therefore, it becomes necessary 

= 7.70 to examine independently the dependency 


LoVer 2.72 


of the magnitude of X on the above- 


mentioned two quantities. . First, the 
0 ¥ 10 1S 20 25 % TiO 

Fig. 4 dependency of. % on grain size was ex- 

perimentally examined under the condition that the content of TiO: is kept constant. 

It has been already known that the initial permeability becomes smaller as the 
dimensions of the grains are diminished, since several experiments upon this subject 

have been reported.(10),(11) In the present study, the specimen No. 1, having the 
largest grain size, was pulverized into the smaller grains successively, and the sus- 

' __ ceptibility of their each step of grain size was measured. The results are shown in 
: Fig. 5 together with those of the same test in an artificial magnetite and of R. 


Chevallier and S. Mathieu’s similar work upon hematite. These three experiments 


may be considered to show nearly the 
same relation between the suscepti- 
bility and grain size. As the dimen- 
sions of the particle are diminished, 
the relative contribution of the various 
energy terms to the total energy are 
changed(7) The wall energy of the 
boundary surfaces between domains is 
a surface energy and the magneto- 
static energy is a volume energy, be- 
ing proportional to r*? and r>_ respect- 
ively, where r is the linedr_ dimen- 
sion of representative grains. 

For the small value of r. then, 


Fig. 5 


the surface energy becomes more dominant than the volume energy. Therefore, 
for the purpose of minimizing the total energy, the total area of the walis must 
decrease at the expense of the increase of volume energy. The magnetization change 
in a weak field is considered to be due to the reversible wall displacement, so that 
it is expected theoretically that the initial permeability becomes lower as the parti- 
cle becomes smaller. The experimental results. shown in Fig. 5 are considered to 
coincide with the above-mentioned theoretical expectation. 

Now, the observed magnitude of X of the fifteen specimens is plotted against 
the corresponding grain size in Fig. 6. 
It is obvious in this figure that the Es 
decreases of X with the decrease of O 
grain size is more rapid than that 


expebted -from: theorestiltgiven By 9°” - : ae 

Biss. SATs fackisrillesusensbothat ee ~ , 
there:ds ark othercelfact: coctrolinee. 2° ' use sosiii 
the magnitude 2 of susceptibility, ; which | eo es ~~ i. 


may be attributed to: the content of ieee Be ingaeel 
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the specimen No. 1. 
The values of %) thus obtained were plotted against the content of TiO. in 
Fig, 7. From these results it may bea summarized that the magnitude of initial 
id susceptibility depends on the content 
a PS of TiO., the former decreasing with the 
° increase of the latter. In the case of the 
a fc) = reversible type specimens in which TiO, 
is assumed to appear chiefly as ilmenite, 
= ° the effect of TiO, upon %, may chiefly re- 
sult from the change in the conteht of 
Fe;0O; in them. On the other hand, in 
the irreversible type specimens in which 
0° the Curie-point decreases with the in- 
crease of the content of TiO. and TiO, 
1 1 pbs ii 
S 


—_" + 
10 1S 20 25 Ti®2% Was assumed to be dissolved into Fe;O, 
Fig. 7 


~) 


and to be present in a complex ferro- 
magnetic mixed-crystal, the effect of TiO. upon %) may be chiefly due to the mag- 
netic property peculiar to that mixed-crystal itself. 

§ 7. Conclusion 

From the experimental facts obtained in the present study, the chemical cons- 
tituion may be considered to play an important role on the magnetic properties, 
especially on the Curie-point and the initial susceptibility of the ferromagnetic 
minerals, though they used to be called simply “magnetite” in the broad sense of 
the word in petrology. The grain size seems to influence the rather secondary 
effect upon the magnitude of susceptibility. 

The reversible and irreversible properties of the thermal change of suscepti- 
bility also seem to be closely connected with the condition under which the original 
rock was formed from an original magna. 

If such problems as mentioned above are resolved completely in the future, it 
could make a contribution not oniy to geomagnetism but also to petrology. 

In concluding this paper, the writer wishes to express his sincere thanks to 
Dr. T. Nagata for his detailed direction throughout this study, and to Dr. H. Kuno 
for his kindness of putting the specimens at the writer’s disposal and of offering 
many valuable suggestions from the petrological point of view. 
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Development and Decay Processes of the Bay Disturbances 
in Geomagnetic Field* 


By Naoshi FUKUSHIMA 
(Geophysical Institute, Tokyo University) 


Abstract 


In this study, the progressive change in the current system for 
the bay disturbance was chiefly examined in four examples with the 
aid of the world-wide magnetogram copies. The average current 
system corresponding to the maximum stage of the bay disturbance, 
which is given in Fig. 3, can be approximated as a result of the 
dynamo-theoretical calculation under the assumption on the distribu- 
tion of the electric conductivity in the upper atmosphere such as 
that the ratio of the conductivity in the auroral zones to that in the 
equatorial region is about ]4 in the dark hemisphere and about 4 in 
the sunlit hemisphere. The development. and decay processes of 
individual bay disturbances are not a simple intensification and 
diminution of the mean current system such as given in Fig. 3. It 
was noticed that the westward auroral zone current in the dark 
hemisphere seems to become most intense nearly at the same time 
when the disturbance force at most part of lower latitude regions 
reaches its maximum magnitude,-and the eastward auroral zone 
current in the sunlit hemisphere does not yet fully develope by that 


time and becomes most intense a few tens minutes later. 


§ 1. Introduction 

The average characteristics of the bay disturbance in geomagnetic field have been 
investigated by many authorities, for example by L. Steiner [1], H. Hatakeyama [2], 
H.C. Siisbee and E.H. Vestine [3], J.M. Princep [4] and others. Some of them intro- 
duced the current system equivalent to the distribution of the disturbance force of 
the bay disturbance with the aid of the world-wide data [2][3]. On the other hand, 
it has been shown by T. Rikitake [5] and further extended by the present writer [6] 
[7] that the current system corresponding to the bay disturbance or to the Sp-field 
can be approximately derived from the dynamo-theory under the appropriate 
assumption on the distribution of the electric conductivity in the earth’s upper 
atmosphere. The investigations of the bay disturbance hitherto, however, concern 
chiefly the average state of the bay disturbance, and so the conclusion obtained 


ought to correspond to its quasi-stationary state. . 


* Contribution from Division of Geomagnetism and Geoelectricity, 
Geophysical Institute, Tokyo University. Series Il, No. 21. 


ial 
\ 
\ 


60 

The polar magnetic disturbances such as Sp and bay disturbance have been 
considered to be caused by the corpuscular stream impinging upon the earth’s upper 
atmosphere after the travel from the sun. Trials for detecting directly the im- 
pinging particles have been made by many observers. The recent results by L. 
Vegard [8], A.B. Meinel [9] and C.W. Gartlein [10] show a strong evidence of the 
presence of protons incoming into the upper atmosphere in high latitude regions at 
the time of intense auroral display. Then, the electric conductivity in the upper 
atmosphere will be fairly affected by such particles at that time, or at the time of 
polar magnetic disturbance, Since the geomagnetic disturbance at high latitude 
regions is considered to be closely related to the ionization of upper atmosphere 
there, knowledges of the characteristics of the polar magnetic disturbances will give 
a clue to research the nature of the corpuscular stream from the sun as well as 
the effect of the impinging particles in the earth’s upper atmosphere. It will be 
significant at the present stage of. the study of the mechanism of the bay disturbance, 
to trace the progressive change in the bay field, or in other words, to know the 
dynamical behaviour of the development and decay of the bay disturbance. 

In the present study, therefore, the whole process of four bay disturbances were 
examined in detail, especially with respect to the characteristics of development and 
decay of their equivalent current systems. 


§ 2. Data used in the analysis 


Examined data were chosen from those occurred in April 1933, during the 
Second International Polar Year Observation, since many magnetograms over the 
world in the equinoctial season could be obtained only for those in this epoch. Dates 
and times of the bay disturbances are 


April 3, 1933, 10—13h G.M.T., 
April 9, 1933, F=12 b°G.M- TE 
April 10, 1933, 14-17h GNLT,, 
April 23, 1933, 5— 8h GM.T. 


The world-wide magnetogram copies at the time of the above four bay disturbances 
were prepared by Prof. T. Nagata at Department of Terrestrial Magnetism, Carnegie 
Institution of Washington, and from U.S. Coast and Geodetic Survey. The number 
of the observatories amounts to 33, 26 of which are in the northern hemisphere and 
4 in the southern one. Besides these data, the tables of hourly values in the Polar 
Year publications could be utilized for reference at 15 observatories. 

The former 33 observatories will be called here the first kind observatories, 
and the latter 15 observatories the second kind ones. The abbreviation and the 
position of the stations and other elements are tabulated in Table I, in which the 
stations written in italic mean those of the second kind. The distribution of the 
observatories is shown in Fig. 1, in which the stations of the first kind are given by 
full circles and those of the second kind by hollow circles, | 
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_ Table I 
Observatory Abbreviation 
p A () x v D D-v 
Thule TH  88°.0 0°.0 76°.5 291°.1 0°.0 -81°.2 -81°.2 
Godhayn GO 79.8 32.5 69.2 306.5 . -175% =5712. <307 
Scoresby Sund SS 75.8 81.8 70.5) ° 383.0 -36.2° "34.5 1.7 
Angmagssalik AN 74.3 52.7 65.6 87:6. A226 4299.8 17 
Sveagruvan SV 73.9 130.7 77.9 168.29 462° 2.9% 247-3 
Calm Bay CB 75) 1533 80.3 52.8 ~ ~92.2 21.1 53.3 
Bjornoja BJ 71.0 124.7 74.5 19:2 37). Saat 36.0 
; Julianehaab Ju 70.8 35.6 2 = 31400 43S 443) 99.6 
5 Fort Rae FR 69.0 290.9 62.8 243.9 24.1 37.5 13.4 
: Tromso TR Gils 1367 69.7 39. -908 37 27.1 
Petsamo PE 64.9 125.8 69.5 S122) - 27.6 5.8 Boe: 
Matotchkin Shar MS 64.8 146.5 73.3 B6A | -22.4°% cay 44.1 
College co 64.5 255.4 649 212.2 27.0 30.4 34 
Sodankyla so b38.< 1205 - “674° 928% 26.7 3.0 29.7 
_ .. Dickson: DI SOS sais e5e) B04) | 128-8. aS 
Kandalkcha KN 62.5 124.2 67.1 32.4 -25.0 6.7 31.7 
Lerwick LE Ss 886s 601. 935885 23.6 -— -136 10.0 
-Dombas DO 62.3 100.0 62.1 GAN 593 6 4 BR ae 
Meanook ME 61.8 301.0 54.6 246.7 172 26.4 9.2 
Sitka SI 60.0 275.4 57.0 224.7 21.4 30.0 8.6 
Eskdalemuir ES 58.5 .. 82:59. - 24553. 3568 <-20A -14.2 6.2 - 
Lové LO 58.1 105.8 EG.t Sea 8.204 96 19 5 
Agincourt _ AG SSO oar B AgT aS So SAN a 96k 07,6 aa eat 2 
Abinger AB 54.0 83.3 512 359.6 -184 -11.5 6.9 
Gross Raum GR 53.3. 105.3 S18. 205 4299.5 .- 251.6 18.5 
Yakutsk YA 51.0 1938 620 129.7 iia 16 San Boe 
Cheltenham CH 50.1 1 3505 -\- 38.7 283°2 2A ATO PEDS 
San Miguel SMe a5 Gan (S-- 2A 38h 113 18 8 Se 
Ebro EB 439° 79.7. 40.8. 0.5  =15.0° . - 9,852 
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§ 3. General description of the geomagnetic condition of the days, on which the bay 
disturbances occurred 
The international daily magnetic character figure C, and Kp-indices on April 
3, 9; 10 and 23, 1933, are given below, where the Kp-indices written in Gothic type 
correspond to the bay disturbances. 


G.M.T, Kp-indices 


Date © | 00-038 03-06" 06-09» 09-12 19-15n 15-18 18-21% 21-24h Sum 
1933 April 3| 0.9 24. ah. 3_ as Va 4s is 24. 2s Awe 
oe oe er i a a BF 135 a i; Tis ie 195 

10,1 0.7 | 3 24 2h, 2 a. Bo 3_ Oe 172 

234 Bi oa Ay Ay. Zo 4 rH 4g Seana 


The mean values of C and daily sum of Kp on international quiet days in 
April 1933 are 0.12 and 7 respectively, while those on international disturbed days 
are 1.28 and 26. The geomagnetic condition on the above four days was not so 
calm, but rather moderately disturbed on the whole. Kp-indices at the time of the 
bay disturbances show a slight increase as will be seen in the above table. 

Some examples of the magnetograms at the time of the bay disturbances are 
reproduced in Fig. 2 (a){d). As will be seen in the figure, the bay disturbance 
observed in high latitude regions is much more complicated in its mode of variation. 
Moreover, the geomagnetic field is not in undisturbed state before and after the 
very time of the bay disturbance in high latitude regions. So that there may be 
a little uncertainty in defining the normal geomagnetic variation which ought to take 
place if the bay disturbance concerned did not occur. It was assumed here that the 
form of the normal geomagnetic variation is a linear change from the value before the 
beginning of a bay to that after the bay. Practically, the said definition of the 
normal variation was applied for 10-13h G.M.T. on April 3, 9-12h G.M.T. on April 
9, 14-17 h G.M.T. on April 10, and 5-8h G.M.T. on April 23. The deviation from the 
above-mentioned normal variation curve was defined as the disturbance force of 


the geomagnetic field due to the bay disturbance. 


$ 4. Character of the mean current system of the four bays at theiy maximum stage 
The disturbance force of the magnetic field at the time of the bay disturbance 
in the low and middle latitude regions seems to increase and decrease simultaneously 
cee the world in a first approximation, though the variation in the disturbance force 
is fairly complicated in high latitude regions. The times at which the disturbance 
force in the low and middle latitude regions reaches its maximum value were deter- 
mined to be 11h 52m _G.M.T. on April 3, 10h 00m on April 9, 15h 00m on April 
10 and 06h 20m on April 23 for respective bay disturbances. The overhead current sys- 
tem ere easanding to the average distribution of disturbance fields at the above-men- 
tioned times is given in Fig. 3, while each current systems at four respective times 
are illustrated later in Figs. 7-10.. The overh2ad current system was derived in the 
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following way. According to T. Rikitake’s result [11], the internal origin part of the 
magnetic disturbance force due to the induced electric current within the earth, is 
nearly a half as much as the external origin part in its magnitude, and the phase 
lag between these two parts is negligibly small. Hence, the overhead current system — 
equivalent to the external origin part of the bay disturbance was so derived that 
the overhead current intensity is responsible for two-thirds of the horizontal com- 
ponent of the disturbance force at respective stations, while the external origin part 
of the vertical component, evaluated to. be twice the observed value, being taken 
into account only qualitatively in the present study. The current system in Fig. 3 

NORTHERN HEMISPHERE IN is in good agreement with Silsbee and 

GEOMAGNETIC COORDINATE Vestine’s one [3], though there is a little 
discrepancy in the phase of the current 
system with respect to the direction to- 
wards the Sun, which amounts to about 
30° in angle. 

The above current system was 
compared with the theoretical one calcu- 
lated by the dynamo-theory [7][12]. The 
calculation was made under the following 
assumptions. The distribution of the 


J 


electric conductivity in the earth’s upper 
atmosphere is assumed to be expressed 
in the geomagnetic coordinate by 


K=ak, in the polar caps 
(9=0°~20° and 160°~180°), 
K=6K) (1+rycos¢) in the auroral zones a 


Fig. 3. Average current system of four bay 


disturbances at their maximum stage. ~ (9=20°~25° and 155°~160°), 4 
. = seactihagecca eee ee en K=K, in the equatorial region | 
direction indicated by arrow. 77 — (@=25°~155" ae 
where K, is a constant value, and the Constante a,b andy are so chosen that =m tl 
3 if anh b>1 and O<r<h., . : eae) wi ¥ 
in the expression of the conductivity in the auroral zones denotes. , 
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In the numerical calculation, only the most predominant term of x~=1 and m=1 was 


9 


sdovted. - Ine ae 2 = ; sh 
adopted, and a! was assumed — ->@ for convenience [6]. Then, the current function 


a 


J derived from the dynamo-theory has a form 
J=Iq +Jo, 
where J, corresponds to S,-field, and Jp to Sp-field attributed to the anomalous large 
conductivity in the auroral zones. J» in the above relation can be separated as 
Jo=Joot+rJ'v, 

where Jp corresponds to the case when y=0, and yJ'v is a perturbation term due to 
the longitudinal inequality of the conductivity in the auroral zones. The details of 
the calculation has already been described in a previous paper [7], together with the 
results of computation. 

The observed currret system of bay disturbance given in Fig. 3 was compared 
with the Jp current systems derived theoretically by varying the magnitude of a, b 
and y. It was concluded then that the distribution of the conductivity in the upper 


atmosphere at the time of the maximum deviation of the bay disturbance force in 


NORTHERN HEMISPHERE IN the low and middle latitude regions can 
GEOMAGNETIC COORDINATE 
Oh be best expressed by 


GN Oe anda 0.0 


in the present model, or in other words, 
the conductivity of the auroral zone is 
about 14 times that of the equatorial 
region near midnight and about 4 times 


near noon, while the conductivity in 
the polar cap is comparable or some- 
what large compared with the equatori- 
al one. In Fig. 4, the current system 
calculated under the assumption that 


a=1.2, b=9 and y=0.6 is given for 
comparison. Since the current intensity 
in the calculation is determined by the 


product of ky! and A), the value of 
Fig. 4 Current system calculated from dynamo- k,{Ky is appropriately chosen, so that 
; theory under the assumption @=1, b=9 ; : 

and y=0.6. the numerical difference of current 
function between the pole and the equator in the both current systems in Fig. 3 and 
Fig. 4 agrees with each other. The calculated current system shows good agreement 
with the observed oue. The distribution of the values of the current function with 


respect to latitude in the noon and midnight meridians is compared in Fig. 5 in the 


- both observed and calculated ones. It will be noticed that the width of the auroral 
‘zone, which is taken to be 5° in latitude in the calculation, may be somewhat wider 


im the observed one. The neglection of the obliquity of the earth’s magnetic axis 
with its rotation axis, and some other simple assumptions on the air motion and. the 
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§ 5. The variation in the : eg ed eee 


disturbance force with 2 
the progress of time 


As already men- 
tioned in §4 and shown Fig. 5 Change in the value of current function with latitude in 
the noon and midnight meridians. The top and the 
bottom of the figure correspond to the average current 
the disturbance force in systems given in Fig. 3 and to the calculated one in 
Fig. 4, respectively. 


in Fig. 2, the variation in 


high latitude regions 
with time is fairly complicated, while it seems that the disturhance forces observed 
at various localities in the low and middle latitude regions increase and decrease 
almost simultaneously with each other. The times at which the horizontal intensity 
and declination of the disturbance force at all stations showed their optimum values 
during respective bay disturbances are plotted in Fig. 6, where the circles show 
those of the horizontal intensity and the squares those of declination. The full circles 
and the black squares correspond to the maximum values of horizontal intensity and 
eastward declination respectively, while the hollow circles and the white squares to the 
minimum values. Large circles or squares and small ones correspond to major and 
minor optimum values. It will be seen in Fig. 6 that the circles and squares do not 
always concentrate around a few definite time-lines, especially for the cases of the 
observatories in high latitudes. This fact shows that the variation in the disturbance 
forces in high latitude regions is so complicated that some of the maximum and 
minimum of magnetic disturbances at different places take place at different times. 
However, there are a few time-lines around which are concentrated most number of 
the optimum symbols, though not all of them. These times are, for example, 

on April 3; 11h 06m, 11h 30m, 11h 52m, 12h 05m, 12h 40m, 

on April 9, 09h 40m, 09h 50m, 10h 00m, 11h 00m, 

on April 10, 14h 30m, 15h 00m, 15h 10m, 15h 30m, 
- on April 23, 05h 48m, 06h 00m, 06h 10m, 05h 20m, 05h 40m, O7h 20m. 
The change in the current system during a course of the bay disturbance will be 


described for respective bays in the followings, and the common characteristics are 
summarized later. 
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(i) The bay disturbance on April 3, 
1933, 10 h-13h G.M.T. 

The current systems at 11h 06m, 
1lh 52m and 12h 05m are illustrated in 
Fig..7 (a)}-(c). 
Figs.8-10, the eletcric current of 1.0 


x10° amp. flows between successive 


In Fig. 7 and also in 


stream lines in the direction indicated 
by the arrow. At 11h 06m, the intense 
westward auroral zone current in the 
dark hemisphere could be clearly seen, 
while no distinct eastward auroral zone 
current appeared in the opposite hemis- 
phere. In the most part of the low and 
middle latitude regions, the disturbance 
force in horizontal intensity was negative 
in its sign, namely the westward zonal 
current was predominant there. At 11h 
52m, the intensity of the westward 
auroral zone current reached its maxi- 
mum. The easward auroral zone cur- 
rent in the sunlit hemisphere was 
formed by that time. In the low and 
middle latitude regions, the disturbance 
force took its maximum value at that 
time in general, and the positive bay 
was observed in almost 40% of the 
whole regions. During 1] h 52m—12h 
05m, the eastward auroral zone current 
scarecely decreased in its intensity, while 
the intensity of the westward one fair- 
ly dimished. In the low and middle 
latitude regions, the area of the region 
in which the positive disturbance force 
of horizontal intensity observed, seemed 
to become a little wider during the in- 
terval of time. This tendency of the 
broadening of the positive bay area, 
which is first noticed by H. Hatakeyama 
(reference [2], p. 24) and later discussed 
by the present writer [14], was clearly 
recognized in this bay disturbance, _ 


Fig. 7 Distribution of the overhead current 
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(ii) The bay disturbance on April 9, 
1933, 09 h- 12h G.M.T. 

The current system at 09h 40m, 
10h 00m and 11h 00m are illustrated 
in Fig. 8 (a){c). At 09h 40m, the west- 
ward auroral zone current in the dark 
hemisphere had developed by that time, 
while the eastward one in the sunlit 
hemisphere was not so definite, just 
similar as in the preceding example. Dur- 
ing 09h 40m-10h 00m, 


auroral zone current became a little 


the westward 


weak in the intensity and the eastward 
current appeared near the auroral zone 
of the sunlit hemisphere, though not so 
intense. The disturbance forces in the 
low and middle latitude regions took 
their maximum values about that time. 
Until 11h 00m, the intensity of the 
current system dimished as a whole, 
while the westward auroral zone current 
became again intense. The broadening 
tendency of the positive bay area, which 
is mentioned with respect to the case 
of the preceding bay, was remarkable in 
this case, too. 

(iii) The bay disturbance on April 10, 
1933, 14h-17h G.M.T. 

The current systems at 15h 00m, 
15h 10m and 15h 30m are illustrated 
in Fig. 9 (a){c). The disturbance force 
reached its maximum value at 15h 00m 
over the whole world except near the 
auroral zone in the sunlit hemisphere, 
where the: disturbance at 15h 10m was 
a little large than that at 15h 00m. In 
other words, the disturbance force began 
to diminish after 15h 00m everywhere 
over the world except only near the 
auroral zone in the sunlit hemisphere. 
The current system at 15h 30m showed 
no distinct vestige of the intense 
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auroral zone current. On the other hand, it 
was also derived that there was no intense 
current system corresponding to the distribu- 
tion of the disturbance force at 14h 30m. 
Therefore, the bay disturbance on April 10 
had comparatively short life, among the four 
bays examined in this study. 

(iv) The bay disturbance on April 23, 1933, 
05h-O8h G.M.T. 

The geomagnetic condition in the high 
latitude regions before and after the bay 
disturbance is not so calm in general. Among 
the four examples of bay disturbance dealt 
with in this study, the one of April 23 was 
in the most good condition, namely the geo- 
magnetic field did not show large perturba- 
tion in high latitude regions before and after. 
the bay disturbance. The mode of the varia- 
tion in the geomagnetic field, however, was 
a rather complicated one. As already shown 
in Fig. 2 (d), the magnetogram in the low and 
middle latitude regions reveals that as if two 
bay disturbances took place successively. 

- The record obtained at Fort Rae showed aiso 
the same character, and the decrease in 
horizontal intensity was anomalously large, 


On the other hand, the records at Godhayn 


and Angmagssalik, which situate south of 
that time, ard the maximum disturbance force 
e at about 07h 20m, at which any 
: ding variation could scarecely be 
Rae and other stations of 
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Fig. 10 Distribution of the overhead cur- 
rent arrows and the equivalent 
current system for the bay dis- 
turbance on April 23, 1933. 
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disturbarce force at the south of Greenlard. 
In the current systems of (a) and (b) of Fig. 
10, we see at once the presence of greatly 
intense auroral zone current in the dark 
hemisphere, while the eastward one in the 
sunlit hemisphere was not yet formed at 05h 
48m and not so distinct even at 06h 20m. 
The current system at 07h 20m showed the 
fadeaway of ihe so-called Sp-current system, 
and there could be seen the presence of some 
disturbance only in the polar region. 

The increase and decrease of the dis- 
turbance force over the world in every 15 
minute intervals during 05h 30m-07h00m 
is schematically shown in Fig. 11 (af). The 
Stations at which the magnitude of the dis- 
turbance increased in 15 minute interval are 
plotted by full circles, and those at which the 
decrease in the magnitude of the disturbance 
force was observed are plotted by hollow 
circles, and the stations at which no significant 
increase or decrease took place are marked 
by the hollow circle with small full circle in it. 
There can be seen a rather systematic dis- 
tribution of these circles. The variation of 
the disturbance force during 05h 45m-—-06h - 
00m given in Fig. 11 (6) chiefly concerns the - 
decrease in the intensity of the westward 
auroral zone current after its maximum at 
05h 48m. The next figure (c) is chiefly due 
to the second increase in the westward auro- 
_ral zone current. Fig. 11 (d), (e) and (/) are 
of the intervals at which the westward 
auroral zone current diminished after the 
second maximum of its intensity at 06h 20m, 
ard the disturbance force became larger at ~ 
the south of Greenland. This series of figure 
illustrates the fact that the disturbance force » 
of the bay disturbance does not: increase or 
decrease simultaneously over the world. 


The is to say, there are cases, in which the | 
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Fig. 1 Distribution of the increasing and decreasing areas in the disturbance force. Stations ie 


at which the disturbance force became larger or smaller with time are plotted by 
@ and © respectively. @ means the station at which no significant increase or 
decrease in the disturbance force was observed. 


disturbance force at some region increase or decrease contrary to the variation-in — 


most part of the whole earth’s surface. The variation in the intensity of the west- 

ward and the eastward auroral zone currents, which were already explained in the 

preceding three examples of the bay disturbances, are of course the case. — 

_ Taking a general survey of the progressive change in the current system of 
the bay Bspurbance during its development Bed ae Process, yne mM : 
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Fig. 12 Change in the total amount of the 
westward electric current in the auro- 
ral zone in the dark hemisphere 
(denoted by W) and the eastward 
one in the sunlit hemisphere (denot- 
ed by E). 
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appropriate assumption on the distribu- 
tion of the electric conductivity in the 


The result 


of the examination of the progressive 


earth’s upper atmosphere. 


change in the current system for the 

respective bay disturbances, however, 

shows that the development and de- 
cay of individual bay disturbance is 

not a simple intensification and dimi- 

nution of the mean current system. 

Summarizing the mode of progressive 

change in the current system of four 

examples of bay disturbances, the 
following facts will be especially 
marked: 

(i) The disturbance force does not 
always. increase or decrease 
simultaneously over the whole 
world, especially in high latitude 
regions. 

(ii) The westward auroral zone cur- 
rent seems to become most 

intense nearly at, the same time 

when the disturbance force at 
most part of lower latitude re- 
maximum 


gions reaches its 


magnitude, and the eastward 
auroral zone current does not 
yet fully develope by that time, 
and becomes most intense a few 
tens minutes later. 

(iii) In the low and middle latitude 
regions, the area in which the 
positive bay is observed, becomes 


wider with time. 


For the purpose of establishing the above-mentioned results, further examination by 
means of many other examples of bay disturbances will be necessary, especially 
with respect to the character (ii). The more detailed examination of the progressive 


future. 
In concluding, the writer wishes to express his hearty thanks to Prof. T. Nagata, 


-change in the current system and consequent conclusions will be reported in the 


who kindly arranged the world-wide data for the present study and gave him kind 
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direction and encouragement throughout the work. Prof. T. Nagata and the writer 
thank to Directors of Department of Terrestrial Magnetism, Carnegie Institution of 
Washington and U.S. Coast and Geodetic Survey, throug nih ey these data 


were prepared. 
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The Effects Of Atmospheric Motion And Dynamo Current 
On The Electron Density Of The Ionosphere 


Ken-ichi MAEDA 


(Electrical Communication Laboratory, Ministry 
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Summary 


The paper deais with the possibility that the ionospheric 
electron density is affeted by the flow of the charged particles (elec- 
trons, negative and positive ions) associated with the atmospheric 
mass motion and dynamo current from dynamo-theoretical aspect. 
The theoretical relationships among the electron density, aimospheric 
motion and dynamo current are deduced and some numerical examples 
as applied to the F2-layer are given. The effect of dynamo current 
is small in comparison to that of atmospheric mass motion, which 
must account for the dynamo current corresponding to the Sq field 
of geomagnetic variation. 

I. Introduction. 

According to the dynamo theory the ionosphere will be the sheet of overhead 
current system, which must be accountable for the diurnal variation field (Sq) of the 
terrestrial magnetism and the theory can be applied to interpret the Sp field with 
some assumptions on the conductivity ef-the ionosphere in the polar regions. 

In the dynamo theory it is considered that the horizontal mass motion of the 
upper atmosphere exists and it produces in the ionosphere urder the influence of 
the vertical field of permanent magnetism of the earth a system of electric current, 
which is responsible for the Sq field. 

By this theory we are led to the idea that the charged particles (electrons, 
negative and positive ions) are made to move with velocities corresponding to the 
mass motion and the overhead current and as the result the ionospheric electron 
density will be affected by such flows.of the charged particles. . me 

Concerning the vertical ionic drift due to a horizontal air current ard its effect 
on the ionosphere D.F. Martyn™ studied the problem on the dynamo theoretical view. 

It is the aim of the present writer to clarify the theoretical relation among 
the electron density of the ionosphere, the horizontal mass motion of the upper 
atmosphere and the overhead current and to show some examples of a deformation 
of the diurnal variation curve of electron density by estimating the order of magni- 
tude of the above-mentioned effects on the basis of informations hitherto known in 
the field of geomagnetism ard ionosphere. 

II, Fundamental equations for the variation of the charged particles, 
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north po. le Fig. 1 shows the spherical 
" coordinates to be used in 
this report and the univer- 
sal time (T), longitude (A), 
colatitude (@) and local time 
(t) are also shown there. 
The following notations are 
used. 


( Greenwich meridian ) 


A=0 


midnight 


u, v: southward and east- 
ward velocities of 


ively, (vertical mo- 
tion is neglected.) 
nu’, vp’: southward and east- 
ward velocities of 
the charged particles 
Fig. 1 Polar coordinates associated with “the 
P:7, 0, t(=2+T), ¢ is reckoned from midnight. current system res- 
pectively, (vertical current in neglected.) 
v=R+h, R: the radius of the earth, 
h: the height of the layer, 
7, M1, Nz: the number density of the electron, the negative and the positive ions 
respectively. 
Suffices 1 and 2 mean the quantity concerning the negative and positive ions ens 
ively and no suffix the electron. 
q: the rate of ion production, _ 
a: the coefficient of recombination, a 
B: the coefficient of attachment, | ; . = 
sr the coefficient of detachmer . a 
‘Then the differential equations coonetng hem and m canbe writen 
as "pee ne sap tn RMT 
_ 700” Sa “ano ‘ 


South pole 


mass motion respect- — 


tres 
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Here the contribution of electron to the overhead current is neglected as the elect- 
ronic conductivity is negligibly small in comparison to the ionic conductivity. Also 


Uy! 


and v,’ are put to be equal to —uw,/ and —»;' respectively and e is the absolute 
value of‘the charge of an electron. 

If we assume the total conductivity of E and F layers or the contribution of 
the both layers as the current sheet to be roughly cqual (this point is discussed in 


the later section), we get the following relations from (4) and (5) for one layer. 


(6); Nyu;' = 1 oy ta oF nv = a: ny oy 
2eZ4 Mtn, rsin Goi’ 2é24 m+n, 00’ 
PSS Se See ee ee 
2ez24 Nitne rsin Gd,’ 2e24 M12 x06 ° 


Next we consider on some approximations concerning the divergence terms in 
(1), (2) ard (8). For instance, in the case of 0 (nv)/r sin 002, the term of space deriva- 
tive of the density xz (that is vdn/r sin 004) is equal to (0n/0t)(v/co), where ¢o (7 sin 
002/0t) is the peripheral velocity of the point under consideration in the ionosphere 
by earth’s rotation. As v/c) can be assumed as much less than unity, the term of 
space derivative of 2 will give a small change to the left hand side of (1) and there- 
fore can he neglected. Also in other cases the terms of space derivatives of n, n, 
and Mz, can be neglected. Therefore the equations (1), (2) and (3) can be written as 
follows, putting s;=2;/(#;+n2) and S:=22/(m,;+n), 


“#C8 ) ae Se =q- an n—Bbu+yny— nf, 
: i dn, 
(9) % = Bu—yn—-Mf+S8iZ, 
(10) oe =q—ann—Mmf —828, 
where a 9 
u v 
(ly. era By aa G02? 
Pesan | 0 oy 1 
(12) &= 262467" lads 00k sing ao 


The ay He (8), (9) and (10): together with (11) and (12) are the fundamental 
equations connecting 2, 2 and w; (density of the charged particles), #, v (mass 
velocity) and_ J (current function). ae 
III. g, the effect of ionic flow peeaeiated with atethett current. 

The magnetic potential (V) and the current function OD) for Sq are given by 


the following equations. 3) 


ee ES (a) a cos{m(A+T)-+e}, 

i coe 8 2n+1 es 
“3 =—— ag odie U1) A+T)+€ 

(14) “Js = 3 Bibs pe = alle RexP'(0) cosm(+ T)+ €3}. 


em and €”. are given in the following table'l, are the eaiica are the mean of ee 
in 1902 and 1905 and the value of 7 was taken as 6,630 km (h= 260 km) for F2 layer. 
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Table I. (e” in 10-5 emu and ¢7, in degree) 


ee see ee: 2 ai 4 5 
a em ent ee | t. et, ee ene, Ct Lee as Om 
‘1 [43.2] 23 | 6.1] 26 |4ti8| 341 | 
ne | +1.8} 241 3.4) 204 | 40.7 197 | 
3 ~~ | |¥o.9) 87 | 2.5] 39|+0.2| 220 
4 +0.4| 2 0.4] 228 +0.9| 84 


values underlined: equinox, others: solstitial inequality. 


_ Using the above valtes J is written as fellows where A+T is replaced by wi. 
(t. is in second and @ is 0.728 x 10sec.) 


Equinox: J (emu)=—9.7 x 10* sin 9 cos 0{cos(wt +26)-+0.68 sin 0 cos (2wt+204) 
e_ +0.35 sin? cos (8w#+39)+0. 10 sin? cos (4w#+228)}, 

“ Solstitial inequality : if (emu) = =—9.7x 10° sin 0{0. 26 cos (wt +23) 

= +0.15 sin @cos (2wé-+ 241)+ 0.07 sin20 cos (Swt-+ 87)+0. 03sin’9 cos (dot +2) 

| ~ £0.11 (5 cos*@—1) cos (ot +341) +0 04 sin 6 (7 cos°6—1) cos (2wt+197) 

oe “£0.01 sin’d (9 cos?d— 1) cos (Bot + 220)-+0.06 sin’ (11 peg cos (dot +84). 


ee oe ee 


For the calculation of g, 2, must be evaluated. z, is defined with the total conduct- 
ivity (A) and the specific conductivity (¢) at the height under consideration, that is 
260km for F2 layer, where o is maximum and the current sheet is assumed to be 
concentrated. 
K 
24=— 
6 


For F2 layer 2, is rosy estimated as 50km. Thus we obtain the flowing 
equation for g. 7 ee See = 


_[Bauinox: ae (ons/em /ses)= —out costo {1 sin of rr antes eat 
: fy TW Vale | ae 


Ce eae feel Spieine cher +106 sind sin Bat) 40.4258 sin? 
| Sosa noua“ _# (ons/emt/see)= ors of Sin 


81 


2.0 OL salle (17) W=3) eP*O)cos 


n=l m=0 


{m(A+ T)+a%\, 


(18) southward Baas 
— ov ae Fae! 7 
“=~ 700 7, ee do 


(Px (0)) cos Ga +T)+a"}, 


(19) eastward eee : 


— oy mm 
oe yatta: Se 


m- m- PO) 
sing 


Joca/ time 


sin {m(A+T)+a7}. 
Fis. 2. Time variation of gXx10 (ions/cm3/sec). 
S: summer solstice, 
E: equinox, cities (#, v) and the vertical 
W: winter solstice. 


The current function (J), velo- 


magnetic field (7) are connected 
by the following equation,” assuming the total conductivity (K) as constant. 


oy 0 /. oy (0 ai 
(20) sin #042 + 00 (sin a0) £ Ser aa >—(vZ) + Su Zsino)\, 
(21) Z=—2Gcos 8. 


Here we assumed the axis of geomagnetism as coincident with that of the earth’s 
rotation. (this point is discussed in the later ee The both sides of (20) are 
calculated as shown below and then k%; and a’? are evaluated. 


(22) ay 0 (sino A) = sin 0S 3 (n useage Re" P"(6) 


sin 6022 = 06 00 n=1m=0 
cos{m(A+ T)+ 7}, 
0 0 : — n kn 
(23) rk} 902) + G5 Zein | = —2GKsind3} 3 ae 
fe (n+2)V (n+m+la—m+1) Pry (0)+(n— (n+ DV (n+ m)(n—m)P™_, (0) 
cos{m(A+ T)+ a7}. 


Mm 


For the harmonics P%,, of J corresponds to Pm of ¥ and ky and az are obtained. 
Next for the harmonics P;, of J corresponds to Pm,, of ¥ and k”,, and a,, are 
obtained. Finally for the harmonics P”,, of J corresponds to P%,, and P34 of ¥ 
and ky; and o™,, must be calculated by using km, and ay,, determined beforehand. 


The equations for k’s and a’s are given as follows. 


eo m “ie ae Je ee a ee oe 
(25) oa a eh Cm VMTN Ge aly = Oi 
(26) _,* cae) (m-+-2)(2m-+5)en» cos{m(A+ T) + €rys } 
gesggeaited oto E 2. cos {m (A+ T)+ary} 
4 hityg (m+2)(M+4V-2MEW 3 cos m(A+-T7) falas) 


2m+7 
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The following table is obtained for k7 and a7. In the calculation 2G was put 
as 34 emu and K as 10-’ emu. 


Table I]. (k”’ in emu and «? in degree) k1=—2.90x10" emu (K=10~7 emu) 


ee eee a 


: - 


MN 
@=77°(Guam) 


ase wave ks a 
7 EERE ald 
we Pare ake 


-fro 
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x29.-409 0=20°(Tromsé) 


(4) 3 6 9 12 15 Simca 24 


local time local time 
Fig. 3-C Time variation of —fX10® (sec—!). Fig. 3-D Time variation of —/106 (sec—). 
S34 summer solstice, scale (1): K=10—‘emu, S: summer solstice, scale (1): AK=10—’emu, 
E: equinox, scale (2): K=0.5xX10-7emu. E: equinox, scale (2): K=0.5x10—‘emu. 
W: winter solstice. W: winter solstice. 


Then we can calculate f (sec—!) by (11), (18) and (19) using the above table. 


/ 
| 


(27). 


es ae 0.66 x 10-8 : 1 es 
Equinox: f(sec™)= K ( sin 0 + sing )oos(ot + 26) + 1.02(3—2cos?6) 


cos(2wt +204) + 0.69sin@(4—3cos?0)cos Bet + 39) +-0.26 sin?0(5—4cos*0)cos (4ot + 228) ; 


—13 y 
Solstitial inequality: f/(sec—')=— Ee cos 0 = “(6 —4cos*0)cos(wt +23) 


+0.39(11—9cos’@)cos(2wt + 241) +0.26sind(19— 16cos*6)cos(3wt+ 87) 
+0.15sin*6(29—25cos?0) cos (4wt + 2) — =a 112cos'0—173cos"0 + 57)cos(wt+ 269) 
—0.27(75cos*#— 108cos70 +-29)cos(2wt + 90) —0.31sin6(132cos*0—181cos?0 

+41) cos (3wé+258)—0.12 sin?6 (637 cos'd —842.cos*0 + 165) cos (4wé + 108)! : 


Calculated examples are shown in Fig.3 for various colatitudes. 


The 


(28) 


calculated equations for # (cm/ sec) and v (cm/sec) are added in the following. 


: et ee ee . 
Equinox: «#(cm/sec)= ay cca {cos (wf +26)+ 1.02 sin@ cos (2wt +204) 


4.0.69 sin?0 cos (wt +39)-+0.26 sin? cos (4wt +228)}, 


» (cm/sec)= 2810 (sin (ot +26)+ 1.02sinOsin Zot +204) + 0.69sin’d sin (Bol 


+39)+0.26 sind sin (4wf+228)}, ; 
Solstitial inequality: a (cm/ sec) = 437% 10" (0.44(2c08°0 —1) cos(wf+23) 
+0.39 sin 6(3cos?0—1)cos(2wt 4241) +0.268in?0 (4 cos’0 —1) cos (8wf+ 87) 
+0.15 sin*6(5cos?0—1)cos (dot +2)-+0.10(28 cos’9—27 cos’0 +3) cos (wt+ 269) 
~ 4.0.27 sin (15cos!0—12cos?0+1)cos(2wt+90)+0.31 sin’@ (22 cos!0—15 cos’0 
+1)cos (3wt +258)-+0.12 sin? (91 cos*8 —54 cos’0+3) cos (at + 108)}, 


v (cm/sec) = AS cos 0{0.44 sin (ot +23)+0.77sinOsin (2ot+241)-+0.79 sin? 6 
sin (Sot +87) +0.60sin*0 sin (4ot-+2)-+0.10(7 c0s?0—3) sin (wt +269)-+0.55sin 0 (3 cos*d 
—1)sin(2wt+90)+0.31 sin?6 (11. cos*# —3) sin (3wf + 258)+0.49 sin*d (13 cos’d_ 
—3)sin (4ot+ 108)}. 
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V. Differential equation of the electron density (F layer) and its approximate solution. 
The diurnal variation of the electron density (#) is given by solving simul- 

taneously (8), (9) and (10), using (16) and (27) for g. and f. But as it is difficult to’ 

obtain a perfect solution for x, we must be satisfied peel rough estimations for some 


limited times of day. 
In order to discuss in some detail, some mathematical approximations must be 


devised. ~ 
n, m and n, are put as follows. 
dn dx 
= = PES bie — =—, 
giekee Ber” di dt 
dn, dy 4 
(29). M=Nioty, Nig= 8 /a epaa pr 7 ae =at* é 4 
dn, dz ; 
y= Nx +2, = fe BET Sew ee are ; 
| 2 was Nx / avr > dt dt” ; 
M, Mio ANd M29 Were defined as the steady state values of electron, negative and 3 


positive ion densities respectively in the case when f and g are put to zero. @ is 
the above-mentioned steady state value of gq. Then. the equations (8), (9) and (10) can 
be written as follows, considering that f is of the order of less than 1/10 of 8 or +. 


(product terms of x, y, z and f are neglected.) oS 


a 6 
ie . « = ase by-the, 
where Baie i mea 
Joie 6 ‘=7, :. 1= aN, ee aoe 
gg ie shigs | %=B, a%=—aMn, Iy=q—q—fr, te. 
Be Be oe hs=4~ ap Pi Se 3 BSE: aompacal 
| The differential es xis given below. 2S aa 


ee each 9% Be peitoteh val - 


ca rN acc oa 


put approximately for daytime not far from noon as follows. 
We assume the Chapman distribution for q and consider q at the level of its 
maximum value. If we denote gn as q at y=0, we get 
ey ee 
sin (0+) ’ 
Q—qo=—I(1+cos wt), (¢ is reckoned from midnight) 


4 =QmnCOS X = qo- 


(34) q =ol sin at, 
q’=w:.I cos wi, 
pe : _cos.d sin 0 
cos 0 sin@ qn= sin (0+6) qo; 


where d is the declination of the sun. Then the equations (32) and (33) can be solved 

and the solution consists of a constant term, two exponential terms and a group of 

sinusoidal function terms, but it is necessary to know at. least three quantities con- 

oe ( an 
r 


cerning # and/or 2 dE: ) at certain times in order to determine the integral 


constant. 


For a rough estimation, neglecting the term of g or g’, the particular solution 
for the following equation 


d*x 


(35) GE — oe AZ + Bx= —Nof' —(7N20 —aNyNo) f 


will give the part ee the effect of f, which should be superposed to the 
usual results (%)) in the case of sole existence of q and q’, that is the general solu- 
tion of 


ad’ x 
dB 


dx 


(36) di 


gt A B= f+ rdq=w): 


For this purpose a particular solution of (35) is presented in the following. 


(B—m’*w’) cos (mot +6m)+moA sin (mot+ 6m) 
(37) a= BIC (mw A) +(B—m'w*)* 
4. 4 Cn (B—m?’o’) sin (mot +0’) —mwA cos (mut+- od ‘m) 
mal (mw A)? +(B—m’*o’)’ 


where C,, and C’, are put as follows, 
4 4 
(38) —MNof! — (7229 — ANN) f =21Cncos (mot +6m)+23C'm sin (mat +0'm). 


(37) can also be used for a particular solution of (36), if C,, and C’, are put 
equal to the amplitudes of cosine and sine components of q/+7(q—qo). 
As n is equal to m+x% +x, we can evaluate a rough variation curve of the 
electron density for the daytime not too far from the noon. 

For the night time variation of*the electron density the foregoing method is 
not applicable. For an approximate treatment we assume the positive ion density (7:2) 
as constant for the time under consideration. Then the differential equation for 2 
is written below. 


— &n 


TE Ons B+ ptamtap ee cami fiatrtam+prpin= Si7g. - 


(39), 
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When we notice that in the coefficients of the above equation the terms containing f 
and f are very small in comparison to P+y+an, and yam, and that the former terms 
are sinusoidal functions of the first order (/* is neglected), the following solution is 
obtained for the first approximation. The term of g serves for the final value of 2, 


but as it is very small, we neglect it for the brevity. 


2 ts a EES 
= Vig 
Miler le ms pert mE \ + Ae i‘. fv , Gtrtan) f+f 


w?+jo(B+r+anz) wo —jo(b+r+anz) f 


where A; and A; are the constants, which are to be determined from the actual 


conditions. In the above calculations the followings are assumed. 
(B+r+am)y>4ran2. 


Unless 7 at night is neglected (here 7 will denote the coefficiect of detachment 
by some non-photoelectric processes), the decay of 2 consists of two trends (expo- 
nential terms) generally, and under the influence of f each trend involves a small 
wavy component, the period and phase of which depend upon those of f and /’ and 
other constants. 

As seen from (40), the first term is the a en decaying component and its 
Wavy component becomes insignificant soon after sunset. The second term pear 
slowly and its wavy crbor does not fade out io rapidly. 


VI. Numerical examples. 

The curves of g and f were shown in Fig. 2 and Fig. 3 in each foregoing sec- 
tion, and as the curves show,.g can be neglected for most cases owing to its small 
order of magnitude. $ ae 

In order to know the variation of 2 we must ‘detennine the values of qo, a B, 
7, 8, and $. As examples for F2 iver the sonsvenee three sets of values were adopted 
for calculation. 


earn =70 ions/cmbfece, Oe aioe B=10-'sec-!, “-y=3-10- sec, 
$1=0.2, : a Se—0:5- oe: ~ ee eee : 
. wollot 6 tba on oD He 2 ae 
a= 70, gg ee = 10, care ae 
Loy Se Ned6) 8% of SS Fal OGY eA NN AN) — ee woe 
sry” ®, ots ave Bs A 3 e Mie =10-", p= 14x 10- a y=07x10" eS 
nd Sz 


fs & Lab BES hh saga ala 
“For the above values pelt eed h . 


eo 


xxi0" 


Fig. 4-A 


Ss: 
ir 


6=54'( Tokyo) 


local time 


Time variation of x x10—4 
(electrons/cm’). 
ease-i.- 6=10-4 y=3X10-*. 
summer Solstice, scale (1): A=10—’emu, 
equinox, scale (2): K=0.5x10—7emu. 


W : winter solstice. 


(2) W) 0=54° (Tokyo) 


local time 


4-C Time variation of x=10-+ 
(electrons/cm’) . 
case 3. 8=14x«10—, 7~=0.7X10-+. 
S: summer solstice, scale (1): K=10—’emu, 
E: equinox, 
W : winter solstice. 


Fig. 


scale (2): K=0.5X10—“emu. 
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e 
(2) 0=54 (Tokyo) 


rallied 
ae 


| 


-4 -2 | t | 
NI E 
-8 -4 i =a 
-l2 -6 
S 
-16 -8 
-20 -10 
3 6 9 12 1$ 18° 2) 
local time 


Fig. 4-B_ Time variation of *«x10—4 


Roe 


(electrons/cm’). 
case 2. ®=y=10-1. 
summer solstice, scale (1): K=10—’emu. 
equinox, scale (2): K=0.510—7emu, 


: winter solstice. 


concave upward with two peaks in the 
daytime and the bottom of the curve 
comes near noon. For comparatively 
small x such as in winter, the single 
peak of ” exists and its time can fall 
in the forenoon, though the usual 
curve without f and g has generally 
its peak in the afternoon. 

For the night time variation the 
wavy components of slowly decaying 
term, that is the second term of (40), 
is significant. The exponent of the 
slowly decaying term (ye7./(8+7+an2)) 
is of the order of (0.3~0.5)x 107 for 
the above three cases and the magni- 
tude of the wavy component is about 
(1.3~1.4)x(2~3)%. If f for the night 


is taken as three times as great, as discussed in (iii) of the next section, the per- 
centage will rise to about 10% of the slowly decaying term. This point is expected 
to account for the night time variation of # in the F layer. 

As for the E layer the similar effects as in the F layer are expected, but the 
order of magnitude of the effects is comparatively less than that in the case of F 
layer. The considerations about this are left for future study. 


— 


VII. Notes and discussions on some associated problems, 
It.must be mentioned here that there are some points to be discussed in the 
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course of this study, which are directly connected to the physical nature of the dynamo 
theoretical phenomena. 

(i) In the foregoing sections the velocities of electrons, negative and positive 
ions in the mass flow of the upper atmosphere were assumed as all equal (w and »v), 
but there may be some questions regarding the velocity of electron flow. If we take 
into consideration the suggestion that the electrons lag behind the general flow of 
of atoms and molecules under the extremely low pressure of the ionosphere, the 
term of f in (8) should be neglected and consequently C and C’ in (32), (83) become 
as follows. 


C=q'+7q —(rtio— amr) f + rane f+(Sir+S,aN) g’ +ramg, 
(41) 


C= +7(4—Go) — (1210-9) f +-ram?S fat. 


And (39) for the night time equation becomes as follows. 

(42) SE (Btrtamtp) Gr + (ram+f (B+am)}n= sire. 
As seen from the comparison of the above equations with (32), (33) and (39), there 
will be no substantial differences in the contents of consideration of a in the fore- 
going sections, even if w and v for the electrons are neglected. 

(ii) Dynamo theory assumes a single layer as the overhead current sheet. It 
seems however to be accepted according to our recent study™ that the F layer can 
serve as a sheet as well as the E layer. If the horizontal mass velocities in the both 
layers are of the same order of magnitude and the same phase, their contributions 
will be roughly equal, as their total conductivities are of the same order of magni- 
tude, although a strict treatment of dividing their parts of contribution is difficult. 

If we adopt the hypothesis by N. Fukushima,” that there exist double sheets 
of current system and the velocities at each sheet is different in magnitude and 
opposite in phase, the velocities “ and v, and consequently f, become greater than 
those treated in the foregoing sections. Also the ionic flow (g) is made greater, and 
may not be simply neglected in some cases. 


(iii) T. Nagata and M. Sugiura® explained the longitudinal inequality of 
solar diurnal variation of geomagnetic field by taking into account the non-coincidence 
of the axes of earth’s rotation and geomagnetism. By this work it can be suggested 
that similar considerations as the above may lead us to a theoretical explanation of 
longitudinal inequality of electron density variation of the ionosphere (F layer). It 
was not difficult to extend the present study on this line, but this idea had to face 
the difficulty that the order of magnitude of g is not sufficient enough to discuss the 
above inequality of the electron density of F layer. In case of the hypothesis of 
double sheets with opposite velocities, a freedom of choice of each velocity is given, 
but it is not allowed to make g too great independently, because f becomes too great — 
for each layer at the same time and its effect on the electron density grows unrea- 
sonably great. 


It was tried to extend the study to the effect of solar disturbed current system 
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for Sp, taking into account the non-coincidence of the both axes, but for the same 
reason as above it will not be reported here. 


(iv) In the equation (20) connecting J, #, v and Z, the total conductivity (4) 
is assumed as constant. But K varies actually in time and space, though its varia- 
tion may not so great as the electron density (7) itself. We put 


(43) K=K,/é, 


“where € is a function of 4 (or #) and 0, and if we neglect the term of space deriva- 


tives of €, we can rewrite the equation (20) as follows. 


(4) indo + ae sin Gp) =e ae (22) +a rai) 


Therefore the velocities « and v hitherto obtained in the foregoing sections are ver- 


tual or effective ones and equivalent to = and ie respectivly. When we denote the 
former z and v as ileg and ve, then follow the equations for oe actual velocities 
and »v, 
: U=Etteg F 
Ww seal yee 
O=E03,.4 
According to our study & is roughly evaluated as 1 at daytime and about 3 at night 
in the middle latitude, with 10-7 emu of K,, which was adopted in our foregoing 
treatments. 
According 3 hia consideration the value of f at see becomes simply 3 times 
as great as that in the foregoing treatments. £3 a. 
(v) It is to be noticed that by the presence of oS and f the. electricity can not 
be strictly neutral in the ionosphere. | The relation i is es by the following 
equation from (8), (9) and (10). : 


“a 


(46) 2 nnn) =F Sime i 


If g and f BEN neg fagibiy ey Pies escent ity is practically neutral een 

= * ss : : z ce K : set: f : ae 
a oe es : 4 ea a 
a References . ta lin fae TPE 


a 
angie Cae 


_ 
irapten Say 


On the Effect of the Earth’s Magnetic Field | 
on the Virtual Height of the Ionosphere . 


By Teruo SATO 


Geophysical Institute, Kyoto University. 


Abstract 


In this paper we calculate the virtual height of the reflection 
point of a radio wave incident in the ionosphere at a frequency of 
0.834 fc, where fc is the critical frequency of the ionosphere, in the 
presence of the earth’s magnetic field. We also calculate the magni- 
tude of the retardation of the wave resulting from the transmission y 
in the lower ionized layer. In these calculations, we take three cases 
that the inclination of the magnetic field vector respectively is 0°, 30° 
and 90°. It is found that the virtual height of the reflection point of 
the wave at a frequency of 0.834fc, which is equal to the actual 
height of the ionosphere in the case of no_magnetic field, does not 
indicate the actual height in the presence of the earth’s magnetic field, 
and the value of the retardation depends on y, which is the ratio of a a 
the gyro-frequency to the wave frequency, and on f;/f2, where f; and _ =e 
f, are the critical frequency of the lower and upper ionized layers . 7 


Ee a a ee 


oe ae 


respectively. 

Further, it is also found that the effect at the earth’s magnetic 
field on both virtual height and retardation is maximum, when the | 
inclination of the earth’s magnetic field vector is 0°, that is, on the ee 
magnetic equator; and minimum when the inclination is 90°, that is ise as 


on the magnetic Boies 


ee iainariaction.. 


91 
using Appleton’s formula. 

In these studies the effect of-the earth’s magnetic field is not considered. If 
this effect is not great, the results of these studies hold. In general, this effect has 
been considered to be small. 

In this paper we calculate the virtual height of the reflection point of the wave 
at a frequency of 0.834 fc, in the case of vertical incidence in the presence of the 
earth’s magnetic field by the numerical operation, assuming a parabolic electron density 
distribution, and for some typical cases of inclinations of the earth’s magnetic field. 
Further, we calculate, for the assumed ratio of the critical frequency of the lower 
ionized layer to that of the upper layer, the magnitude of retardation of the incident 


wave in the F; and F, layers resulting from the transmission of the wave in the 
lower ionized iayers. 


2. Calculation of Virtual Height in the Presence of the Earth’s Magnetic Field. 


The virtual height of a vertical incident wave in the presence of the earth’s 
magnetic field is clearly given by 


ev 


sot a -1{(19 pitt. Bete) te (1) 


where f is the frequecy of the incident wave. The distance between the earth and 
the base of the ionized layer is omitted for simplicity. 

In order to calculate this integrai, it is more convenient to use the relation 
between the refractive index and the electron density and also the electron density 
and the height, than to use the refractive index itself which contains earth’s magnetic 
field in a complicated form. From the definition of the height, this height is equal 
to the half-thickness of the ionized layer. At first, the relation between the electron 


density and the refractive index is given by Booker as follows ; 


8+ Ax*+ Bxtr=0, (2) 
A=—u{1+XxC—q)}, (3) 

=(C?—¢’) (20 —y") + (Sin + Sz + gyn)? +? (C?—@?)}, (4) 
P=—u{w—y}(C—@_y, (5) 


where x=47e’N/mp’, 
p=2nf, N=electron density, .c, m=change and mass of electron. 
u=1—iz=1, 
z=v/p, é 
S,, Sz, C=components of a unit vector in the direction of incidence along the 
rectangular axes 1, 2, 3. In the vertical incidence, Si=S:=0, Cals 
Pe upward vertical component of phase-propagation vector. 
This equals to the refractive index in the vertical incidence. 
y=fn/f,= ratio of gyro-frequency to wave frequency. ~ 
Its 2, Ys.= Component of the earth’s magnetic field vector along the axes 1, 2,3. 
If we take the axis 3 in the direction of propapation of the wave, the direc- 
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tions of the axes 1 and 2 are arbitrary. Because the components of the earth’s 
magnetic field vector is contained in (Syi+Sx2+qvs) in (4), and this is the scalar 
product of the phase propagation vector and the earth’s magnetic field vector. In 
other words, the effect of the earth s magnetic field depends only on an angle between 
the direction of the phase-propagation and the earth’s magnetic field vector. In the 
vertical incidence, $,=S:=0, q=, and the effect of the earth’s magnetic field depends 
on an inclination of the field vector. As the component ys is contained in the form 
of ys’, results of calculations in the north hemisphere are equal to results of calcu- 
lations in the south hemisphere. In this. study we take three cases in which the 
inclination of the earth’s magnetic field vector is 0°, 30° and 90°. The inclination 0° 
and 90° are the value on the magnetic equator and magnectic pole. In these places, 
the solution of (2) shows a kind of a singular solution. Therefore, it seems to ba 


appropriate to start from the calculation in the inclination of 30°. 


A, Case of Inclination 30°. 
Equation (2) has three roots. One of them refers to the ordinary wave and 


other two roots to the extraordinary waves.. One of the two roots corresponding to | 


the extraordinary waves is. of no practical interest, because its does not show x=0 
for the value w=1, but the other chracterises that the wave is reflected from the 
rn x=1-—y. For a root corresponding to the ordinary and extraordinary wave 
as if and tor is given by 


=| OB. or , Ox Ox? Ox 


POE ct LOT ee Ole denials OF eae Ee , 
of | PAs OB Or s? (6) 
(On _ 2 ont a )(6— —4y?— 2p + p2y?) #4432 (223 +628 
Tap Ha + by? BH) I (Binge 
+ Matte 


jerKor a a ame value of y, (7) is ceateied 39 substitution of the value of x | 
and y« which is related by equation (2). 
rt x we substitute (7) into (1), then ae is. 
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than that of an ordinary wave for the same layer. Curves also indicate that as y 
increases, ee oF increases for the same x. This implies that the effect is great as 
the wave frequency is small. Further, it is known that as the electron density 
approaches the maximum value, this value rapidly increases. Therefore, the effect 
of the earth’s magnetic field is great in the neighbourhood of the refiction point and 
not in the region of the smail electron density. 


As the electron distribution, we take a parabolic form given by the expression 


s=100 {1—(1—x)"}, (8) 


where s is the height in km (that is, the distance from the base of the layer), and 
x=4re"N/mp’. 
This expression shows that at x=1 the height is 100km and at x=0 the height 


is 0. As above mentioned, the sending and receiving station of the wave are situated 
tee 0¢ 


The electron density at a height of the reflection point of the wave at a fre- 
quency of 0.834 fc in an ionized layer, is equal to 0.695 N,, where N, is the maximum 
electron density of the layer. Therefore, the distribution of the electron density to 
be applied to the calculation is the part below this level. Thus we can integrate 
numerically (1) by use of (7) and (8). These results of calculations for the ordinary 
and extraordinary waves are shown in Fig. 1. hn and h, show the virtual height of 


the wave at a frequency of 0.834 fc in the 


1 
— ORDINARY WAVE 
--- EXTRA-ORDINARY WAVE 


presence of the earth’s magnetic field and 
that in no magnetic field respectively. 
(4im—h,) is the value to be corrected in 


the calculation in no magnetic field. As 


— 
a 


(hw ha) IN KM 
Ns 


we presume,these curves explain that the g ag 

virtual height of the wave at a frequen- Fig. 2 The difference of the height of the 

cy of 0.834fc in the presence of the reflection point of the wave incident 
; 5 in the ionosphere in the presence of 

earth’s magnetic field, is greater than thevcanicaduci nce neler 

that in the case of no magnetic field. magnetic field. Am shows the height 

For example, in the case of y=1/4 for an . in the presence of magnetic field and 


ha in no magnetic field. 
ordinary wave, the actual height should : e 


be smaller by 20.5km than that calculated in no magnetic field. Tieratones when we 
calculate the actual height of the maximum electron density, we must take the 
virtual height of the wave at a freqnency less than 0.834 fe. Curves also show that 
the quantity (4m—f,) is greater for small value of 1/y than for great value, and 
greater for the extraordinary wave than for the ordinary wave. It is well known 
that the difference between the critical frequency of the ordinary wave and that of 
the extraordinary wave, is about 0.6 Mc/sec. in practical observation for middle lati- 
tude. Thus (n—ha) for the ordinary wave, for example, at a 4 Mc/sec. cofresponds 
to (4n—h,) for tne extraordinary wave at a frequency of about 4.6 Mc/sec., and (hn 
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—h,) from the two waves is equal. It is sufficient to calculate for only one kind of 


wave. 
These results are proportional to the semi-thickness of the ionized layer. 


B. Case of Inclination 90°. 

Equation (2) has three roots. In this case one of them is 1. This implies that 
the refractive index is independent of the electron density. Then this root cannot 
be taken in our study. Two other roots are (l—y’)(1+y), (l—y?)(1—y). Waves 
characterized by these roots are the wave reflected from the point of x=1+y and 
x=1-y respectively. These refer to the extraordinary wave. roe for two waves 
is given by 

On _ 29? 1=v*) +0. — n°) 2x (8-1? —2y*p*)— 4a" (3—2n*) +68 (9), 
2 (x—1)[*+—-Q—y") (1—2)] 


Curves w+ fot OF for two waves, after a similar operation with that in A, are given 
= in Fig. 3 for y=1/4 and 1/10. Full lines 

show the wave reflected from x=1+y, and 
dotted lines show the wave reflected from 
_x=1-y. It is easily clear that the form 
of curves.are similar with those in Fig. 1, 
and “+ ree for the wave reflected from 
x=1-—y is greater than that for the wave 
reflected from x=1+y. Comparing the 
dotted line for y=1/4 and that for 1/10 in 


wt Fe is greater for lines in Fig. 1 than 


of the earth’s magnetic field becomes 


Fig. 3 The Beles between cere and x 
for the ae ihe poidt greater as the inclination decreases. This 


- x=1-y and PERF for y=1/4 and tendency is indicated over all x, and — 
1/10 in the case of inclination 90°. “ specially remarkable near x=1. <i 
ceTntecratel height.ds aiven.in Fieke. fee aston Se ea 


GA sorter toys ene: cetted ne etowstttat for the: 
asl-y. (lem Ite) for the two waves are ee ge 


Fig. 1 and 3, it is found that the value — 


that in Fig. 3. This implies that the effect 


re wr oe 


wk ere ea 
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yy PAE Gl CD 
x= is ear we )+h é (10) 


For the wave characterized by this root, re is given by 


pole =, HX tO — 2") (yy —p*) x+y" (1— py") 


oF ire (11) 


On a et 
Curves for urls are shown in Fig. 5 for y=1/4 and 1/10. The wave characterized 
by (10) is that reflected from the point of maximum electron density defined by x= 
1—y. Thus comparing the curves in Fig. 5 with those in Fig. 1 and 3, erty in 


Fig. 5 is the greatest in three cases for the same electron density and the same 


Fig. 6 (ha —hm) for wave reflected from the 
point *=1—y in the case of inclina- 
tion 0°. 


value of y. Especially, the difference is 
great near x=1. 


; Integrated height is given in Fig. 6. 
Fig. 5 The relation between Fas ge and x (h,—h,) for the same 1/y, is the greatest 
for y=1/4 and 1/10 in the case of in three cases, _ ; 
inclination 0°. : ; 
After all, the effect of the earth’s 
magnetic field is maximum in the case that the inclination of the magnetic field vector 


is 0°, that is, on the magnetic equator and minumum on the magnetic pole. 


3. Calculation of Retardation. 


By use of the above mentioned numerical we will consider the magnitude of 
the retardation of the wave in the lower ionized layer. It is pointed out that in the 
case of vertical incidence of the wave, for example, in the F, layer, the wave is 
affected by the F; and E layers. In the calculation of this kind it is necessary to 
know the electron density and the thickness of each layer and the form of the valley’ 
between two layers. The form of the valley between F; and F; and that between 
F; and E are not distant. These are not detected by the wave observation. Man- 
ning“ assumed this form in his work. Therefore, we must, too, assume the form of 
the valley in the calculation. We classify this form into two cases, that is; 

(a) the case where the electron density over the valley is the same as the | 
‘maximum electron density of the lower ionized layer, and 

(b) the case where is no another electron distribution particularly defined 
between the two layers, which is superposed on the electron distribution of two largers. 


a aD 
\ 
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A. Case of Inclination 30°. 


(i) Case of the crititical frequency of the F» layer. 
In this case, as the electron density in the E layer is too small compared with 
that in F, layer, we can neglect the effect of the retardadation in the E layer. The 


electron density of the F, layer is not much great but two times of that of the E 


layer. Concerning the case (b), the effect is small compared with the case (a), so 


that the calculation in the case (b) is omitted. The electron distribution of F; and 


F, layer is assumed to be Chapman’s and parabolic respectively. The scale height 


and the height of the maximum electron deusity in the F; and F; layer are 30km, 
70km and 220km, 320km respectively. Results of these numerical calculations only for 
the ordinary wave for y=1/6, 1/8, and 1/10 are shown in Fig. 7, as the magnitude 


of the retardation for the ordidary and extraordinary waves is almost equal. 


These curves indicate that the re- 
tardation of the wave depeuds on the 
ratio of the critical frequency of F; layer to 
that of F., layer. 
the retardation increases. 


As the ratio increases 
Therefore, if 
the ratio of the electron density of the 
F, layer to that of F, layer is small, the 
effect of the retardation can be neglected. 
From the wave observation, it is found 
that the maximum electron density of the 
F, and F, layer in day-time is about 2~ 
3x10°/c.c. and 1~2x10°%/c.c.. Then the 
effect of the retardation on the virtual 
height can be neglected in practice. 
(ii) Case of the critical frequency of the 
F, layer. 

In this case the scale height of the 
E layer is taken to be 10km, and as the 
difference between the electron density in 
two layers is small, it appears that the 


QUANTITY OF RETARDATION IN KM 
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Fig. 7 The magnitude of the retardation of 
the wave at the critical frequency of 
the F: layer in the case of inclination 
DUP s 


retardation is considerable. From the wave observation, it is found that the maxi- ~ 
mum electron density of the E layer in day-time is 1~1.5x10°/c.c.. . 


case (a) 


Results are shown in Fig. 8 for y=1/3 and 1/4 by full lines. It is found that | 
as the ratio of the critical frequency of the E layer to that of the F; layer increases, 
the retardation in E layer rapidly increases as in the case (i), and reaches to 18.5km — 


for 0.7. 
case (b) 


Results are shown in Fig. 8 for y=1/3 and 1/4 by dotted lines. It is clear that 
the effects are smaller than those in case (a). In the practical observations the | 


possibility may be greater in the case (b) 
than in the case (a). Thus when we are 


concerned with a measurement of the 
height of the maximum electron density 
of F, layer, it is necessary to correct. 

B. Case of Inclination 90°. 

(i) Case of the critical frequency 
of the F, layer. 

As the quantity of retardation for 
the two waves are equal, operation for 
the wave reflected from the point x=1—y, 
is carried out for y=1/6, 1/8 and 1/10. 
The quantity of the retardation in Fig. 9 
is small for all y compared with Fig. 7. 
When the ratio of the critical frequency 
of F,; layer to that of the F. layer is small, 


the retardation does not change for the 
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QUANTITY OF RETARDATION IN KM 


(ov) 94) 


Fig. 8 Quantity of the retardation of the 
wave at a critical frequency of the Fy 
layer. Full lines show the case (a) 
and dotted lines the case (b) in the 
case of inclination 30°. 


increase of the inclination of the earth’s magnetic field, and it becomes greater for 


the smaller inclination when the rotio increases 


(ii) Case of the critical frequency of the F; layer. 


Curves for retardation for two cases (a) and (b) are shown in Fig. 10. Full 


and dotted lines show case (a) and case (b) respectively. 


When ratio of the critical 


frequency of the E layer to that of the F; layer is great, the decrease of the quantity 
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Fig. 9 Quantity of the retardation of the 
; - wave at a critical frequency of the 
F, layer in the case of inclination 90°, 


~ QUANTITY OF RETARDATION IN KM 
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_ 
kEALF, 
Fig. 10 Quantity of the retardation of the 

wave at a critical frequency of the 
F, layer in the case of inclination 
90°. : 
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of the retardation, comparing with Fig. 8, is remarkable. This decrease is mainly 
due to the decrease of retardation above the maximum electron density of the E 


layer. 
8 


C. Case of Incliuation 0°. 

(i) Case of the critical frequency 
of the F», layer. ; 

In this case, the quantity of the 
retardation is great compared with Aaud ~ 
B. Fig. 11 shows these quantities. Shapes — 
of curves are similar with (i) in A aud B. 

(ii) Case of the critical frequency 
of the F; iayer. 

Retardation of the wave in the E 
layer is greatest of the three cases. Even 
in the case (b) the retardation reaches, in 

Fig. 12, about 10km for y=1/3 and 1/4 
0! when the ratio of the critical frequency © 
of the E layer to that of F; layer is 0.7. 


% 
: 
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Fig. 11 Quantity of the retardation of the 
wave at a critical frequency of the Therefore, the correction is necessary in 

F. layer in the case of inclination 0°. a measurement of the virtual height of 4 


the F, layer near the magnetic equator. When the ratio is below 0.4, the retardation | 


is small. 


4. Conclusion. 

It is found that the effect of the earth’s magnetic field and the retardation of 
the wave, by which the virtual height of the reflection point of the wave is woe 
increase as the frequency of the wave 
incident in the ionosphere decreases and 
also the retardation of the wave increases 
as the ratio of the critical frequency of 
the lower ionized layer to that of the 
upper ionized layer incresaesy Further, it 
is pointed out that the effect of the earth's 
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effect of the earth’s magnetic field in the layer itself, and the other for the retarda- 
tion of the wave resulting from’ the trausmission in a lower ionized layer. These 
corrections are to be taken in a measurement conceruiug the F, layer, but one of 
the two, the retardation in the lower ionized layer, can be neglected concerning the 
F, layer. 

In conclusion, the writer wishes to express his hearty thanks to Prof. M. 
Hasegawa, Dr. M. Ota and Mr. Hirono for their kind advice, aud to Dr. K. Maeda 
for helpful suggestions during the course of this study. 
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On the Relation between the Cosmic Ray Intensity 
and the Geomagnetic Storm 


By Kazuo NAGASHIMA 
Geophysical Institute, Kyoto University 


Abstract 


Many investigators made various attempts which were calculat- 
ed to explain the decrease of cosmic ray intensity during a magnetic 
storm in the maguetic field of the ring current, suggested by S. Chap- 
man.“ Unfortunately their results seemed to contain many qualita- 
tive faults. 

In this paper, I do not try to search for the cause of the 
decrease in the magnetic field, but in the electric field, connected with 
the magnetic storm. From such a standpoint, the shape of the west- 
ward current system which is appropriate to account for the cosmic 
ray phenomena during the maguetic storm, is discussed. Qualitative- 
ly, the results, obtained, seem to be adequate enough to explain quite 
reasonably the observed phenomena of cosmic rays, although, quanti- 
tatively, there still remain somewhat questionable points. 

In this paper, I also discuss relative merits between the mag- 
netic field hypothesis and this electric one. 


§ 1. Introduction 
The decrease of cosmic ray intensity during the magnetic storm is a world-wide 
phenomenon, which was first discovered a S.E. Forbush and thereafter ‘confirmed _ 
by many others. _ — al 
Its principal characteristics are as follows: ieee a 


as (A) Decrease of cosmic ray intensity —When a magnetic storm breaks 
ee ee arto 
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component during some magnetic storms, the intensities of cosmic rays do not vary 
at all. 

(E) Increase of the amplitude of the diurnal variation of cosmic rays during 
the magnetic storm :—Y. Sekido pointed out that the amplitude of cosmic rays 
increased during the magnetic storm and the curve of this increase had no similarity 
with that of the earth’s magnetic horizontal component or that of the cosmic ray intensity. 

(F) Advancement of the phase of the diurnal variation of cosmic rays during 
the magnetic storm :—At the same time, Y. Sekido showed that the phase of the 
diurnal variation of cosmic rays advanced during the magnetic storm and the curve 
of this advancement had some similarity with that of the earth’s magnetic horizontal 
component or that of the cosmic ray intensity. 

S. Chapman concluded that the axially symmetric part of the magnetic storm 
field was due to the electric current which flowed westward above the earth, a few 
earth radii apart. At the same time, he™ suggested that, in the inner region of such 
current system, the earth’s magnetic field was weakened, while it was strengthened 
in the outer region. This can be regarded as the increase of the earth’s magnetic 
dipole moment in the outer region, so that the cosmic ray intensity was decreased 
during the magnetic storm, by the weil known theory of Lemaitre and Vallarta.” 

Accepting this assumption as true, let us examine in what follows whether this 
magnetic hypothesis is suitable in explaining reasonably the above mentioned pheno- 
mena. 

(i) T.H. Johnson, analyzing this storm effect, concluded that the magnetic 
storm effect of cosmic rays was of the order three hundred times too large to be 
accounted for as an increase of the magnetic dipole moment according to Chapman’s 
hypothesis. 

(ii) As for experimental evidence (B), there is nothing to be said about the 
propriety of this theory, as the motion of a charged particle is not clear at present 
in the earth’s magnetic dipole field and the magnetic field produced by the west-ward 
current system. 

(iii) We cannot explain the experimental evidence (C) because the increase of 
the earth’s magnetic dipole moment seems to give no influence in a region well above 
the “knee” of the latitude cut-off which is supposed to be due to the presence of 
the sun’s magnetic dipole field, suggested by L. Janossy. But, as there is at present 
no agreement about Janossy’s suggestion, it may be unsuitable to judge from this 
phenomenon that this hypothesis encounters with a serious difficulty. 

(iv) Experimental evidence (D) seems to be explained by this hypothesis as 
follows: The larger the distance to the current system (for example Ring Current) 
from the earth’s centre, the larger the magnetic dipole moment of the earth. Then 
the decrease of cosmic ray intensity becomes larger. That is to say, this phenomenon 
depends upon the dimension of the current system. . . 

(v) Experimental evidence (F) seems to be explainable by this hypothesis. 


But it is difficuit to acount for experimental evidence (E), because it is low energy 
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region of cosmic rays, that is prevented to arrive at the earth by the increase of the 
earth’s magnetic dipole moment, and therefore the amplitude of the diurnal variation 
of cosmic ray during the magnetic storm must be decreased. 

Some criticism against this Chapman’s hypothesis is brought forward as follows: 
As pointed out by T.H. Johnson,“ and J. Clay and E.M. Bruin,” the magnetic field 
of the region, in which cosmic rays suffer extremely the influence of the earth’s 
magnetic field, is weakened by the ring current of the magnetic storm, and therefore 
there must result an increase of cosmic ray intensity. This idea is inferred from 
the Bruche’s experiment“ of the ring current modei, which indicates the southward 
shifting of auroral zone. 

T. Nagata, treating mathematically the problem of charged particles in the 
magnetic fields of the earth’s dipole and a ring current, confirmed the Bruche’s ex- 
periment. 

From these considerations, we cannot understand why the cosmic ray intensity 
is decreased by such a ring current which causes the southward shifting of auroral 
zone, because cosmic rays are only different in energy from auroral particles, and 
their motions in such a magnetic field must have some similarity. 

Recently, all attempts by many investigators,¢” following Chapman’s suggestion, 
resulted in failures. 

From the above discussions, it is not appropriate to account for the decrease 
of cosmic ray intensity by this hypothesis. But, supposing from their similarity in 
phase between the decrease of magnetic horizontal component and that of cosmic 
ray intensity, we cannot suppose, these phenomena are produced by quite different 
causes. ‘ 

On the other hand, H. Alfvén®? adovocated a hypothesis that the intensity of 
cosmic rayS was decrease during the magnetic storm by the electric field produced 
in the stream which was emitted from the sun. In the next, I will state the outline 
of this hypothesis and attempt to examine the validity of it. 

A neutral ionized stream emitted from the sun, passing through the solar 
magnetic field, positive and negative charges are accumulated at the boundary of 
this stream by the influence of the magnetic field and then an electric field is pro- 
duced in the stream, as shown in Fig. 1. As pointed out by H. Alfvén, the charges 
accumulated on the sides of the stream tend to produce a current, for example along 
the line AEB. However, the inductance of this loop is so big that its time constant 
t is very long, possibly of the order of years. When the time during which the 
stream has been on is short compared to rt, there is no appreciable electric field 
outside the stream, because the field from the charges on the sides of the stream 
is neutralized by an induced field. Then, before the magnetic storm, for example 
when the earth locates at a point B in Fig. 1, the cosmic rays, which arrive on the 
earth passing through the stream, gain energies corresponding to the potential dif- 
ference of the stream. While the cosmic rays, not passing through it, do not change 
in energies. Therefore the intensity of cosmic rays becomes higher than the normal 


103 


state on the average. After the storm, when the earth locates at a point A in Fig. 
1, the cosmic rays, passing through the stream lose their energies, and not passing 
through it, they make no change in energies. So that the intensity of cosmic rays 
becomes lower than the normal state, on the average. He stated also that the 
electric field of the order of about 10° volts was necessary in order to produce 
a variation of 1 per cent of the intensity. This is a rough review of his hypothesis. 
According to this hydothesis, the intensity of cosmic rays must increase before 
and decrease after the storm, and the quantity of the increase must be almost 
the same as that of the decrease. However, in the cases of the cosmic ray storms, 
commonly observed, the increase does not occur so that this is a difficult point of 
this hypothesis. Moreover, the data of cosmic ray storm which he cited in order to 
explain his hypothesis, is that of the extraordinary increase of cosmic ray intensity 
which seldom occurs. Even in this case, the cause of the increase is not the same 
as that of the decrease, because the increases produces extremely the latitude effect 
at various observatories, while the decrease does not. So that, although his hypo- 
thesis is very suggestive, it must encounter with a serious difficulty. Moreover, he 
did not calculate accurately the influence of the electric field upon the cosmic ray 
intensity. 

From the above various considerations, it is possible to say that the cosmic 
ray storm is produced by some cause growing with the west-ward current system, 
if this cause will explain reasonably 
the above mentioned phenomena. The 
magnetic field of the current system 
may become this cause if Chapman’s 
hypothesis will be improved, but I will 
attempt in what follows at explaining 
these phenomena by the electric field 
(which differs from that of the stream). 

Prior to this, we must examine 
the influence of an electric field upon 
the cosmic ray intensity. 


§ 2. On the Motion of Charged 
Particles in the Electro-Magnetic Field 


In a magnetic field, from the 


D 


\ extended Liouville’s theorem, the 
EARTH intensity of cosmic rays in an allowed 
Fig. 1 Polarization of the neutral ionized direction, defined as the number of 


stream by the solar magnetic field. particles of given energy crossing unit 


7 ali angle per unit time, is the same as it is at infinity and therefore, if the dis- 


> 


4 


tribution of cosmic particles at infinity, is assumed to be isotropic, the intensity is 


precisely the same in all allowed directions. So that we can easily obtain a direc- 
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tional intensity of cosmic rays at any point on the earth by only knowing the inten 
sity spectrum of primary cosmic rays and the minimum energy at that point (cut-off 
energy). 

On the other hand, there is a different appearance in an electro-magnetic field. 
The intensity of cosmic rays in space“ is 


3 
I=+«D, (2.1) 


where m and P are mass and momentum of cosmic ray respectively, and not cons- 


tants, and 


oN 
OV 


ON is the number of particles in dV. From the extended Liouville’s theorem, D is 


D= with dV=0dx-dy-dz-0P,:0Py-6P:. (2.2) 


constant in an electro-magnetic field,“ but P and m are varied from point to point 


by the influence of the electric field, then the intensity changes as follows, 


4I=D-(> Ais -am ), (2.3) 
m m 
From Eqs. (2.1) and (2.3) 

Al_, AP Am __ 3W-4W AW 

Wey ee m  W-(mcece . W’ (2.4) 


where m, and W are rest mass and total energy of the particle respectively and ¢ 
is- the velocity of light. With the relation W=E+m)c?, Eq. (2.4) becomes 


4I=L(E)-4E-T, (2.5) 
where 
9 2 (moe? oe?) (C28 
L(E)= re eas ee — er Ps : (2.6) 


Then, the intensity of cosmic particles changes from J to J—4J, when its kinetic 
energy E decreases by an amount 4E. This is most important fact for the explanation 
of the decrease of intensity, as to be mentioned in the next section. Thus, taking 
this effect into consideration, we can obtain a directional intensity at a given point 
on the earth by getting the cut-off energy and intensity spectrum of primary cosmic 
rays, as well as in the case of the magnetic field. 

The equation of a charged particle in the electro-magnetic field is 


SG (mV) =e +2 VHi. (2.7) 


As the current system produced during the magnetic storm seems to be symmetry 
about the axis of the earth’s magnetic dipole moment, so the electro-magnetic field 
produced by this system is also symmetry about this axis. Moreover, let us assume 
this field to be static.2» ‘ 

It is most desirable to examine the orbit of a charged particle in such a field 
and the earth’s magnetic dipole field, and how the cut-off energy varies at a given 
point on the earth. Unfortunately, Eq. (2.7) cannot be solved analytically and so we 


are gratified with the first approximotion as follows: Some variation of cut-off 
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energ, caused by the electro-magnetic field of the current system, seems to be too 
small to strongly influence the variation of cosmic ray intensity on the earth’s sur- 
face. So we will neglect the influence of the variation of the cut-off energy and 
use that of the earth’s magnetic dipole field.C 


$3 Variation of the Intensity Spectrum and Decrease of the Cosmic Ray 
Intensity by the Electric Field 
Let us suppose that primary cosmic rays are all positive protons, whose 
differential intensity spectrum is J(E)dE, where E is the energy of cosmic rays. From 
the investigations by many workers, next two types are adovocated as the spectrum. 


T(E) dE=J-E~-dE, (3.1) 
I(E) dE=J-e-** -dE. (3.2)8»@2) 


and 


In the present state, as it is impossible to decide which is the true representation, 
we might as well adopt both.®® 

In a static electro-magnetic field, mentioned in Section 2, the variation of the 
energy of cosmic rays at any point in space is determined by the scalar potential 
V alone, and does not depend upon their paths. Let us assume a potential at any 
point on the earth to be Vy,, then, arriving at this point from infinity, cosmic ray 
suffers an energy loss 4E, where 

AE=eV,. (3.3) 
So, the primary energy EF changes to E—4JE. 

Thus, taking Eq. (2.5) into consideration, the intensity spectrum I(£) dE varies 
to (c.f. Fig. 2) : : 
I(E+4E):{1—L(E+4E)- 4E}-dE. (3.4) 
If we neglect the variation of the cut-off energy of cosmic rays at any point of the 
earth, as mentioned in the previous section, the decrease of cosmic ray intensity in 
a given direction by this variation of spectrum is easily obtained at the top of the 


atmosphere, i.e., 
(*, (E+ 4E)-{1—L (E+ 4E)-4E} dE 


x 100, (3.5) 
\p,J(B)dE 


yo %)=| 1— 


where E, is the cut-off energy in a given direction at geomagnetic latitude 2°. 
But on the earth’s surface, how the intensity of cosmic rays varies is not 
simply determinable by the atmospheric effects. So I will estimate this order in the 


next. 
[A] In the case of J(E)=J-E~*-dE 
(i) From Eq. (3.5), at the top of the atmosphere, the one of intensity 


becomes 
aie = =r) | 
ye7é)=100-| 1 (ert 


4E Ey [2-2 (4m) 
+100 (a Ss )( ere abi Brae 
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2 z 2, Sos 
+3, aD ) $e i ) Pasi | (3.6) 

(ii) The decrease of intensity on the ground:—In the case of the variation of 
the earth’s magnetic field, there yields only the variation of the cut-off energy, so 
we can make the order estimation of the decrease from the latitude effect of the 
field-sensitive rays at sea-level. 

But in the present case, the situation is slightly different from the former, due 
to variation of the differential intensity spectrum, as in Fig. 2. Let us make two 
assumptions as follows: 

(a) From the latitude effect, when field-sensitive primary cosmic rays with an 
interval of 1 Bev. arrive on the ground, they correspond to 0.7 per cent of the total 
intensity on the ground. Let us expand this experimental evidence up to the region 
IE) of non-field-sensitive rays. 

; J(E+AE)(1- L(E+ AE)-AE) (b) As a result of this expansion, 
ag the primary cosmic rays have an effect of 
he 


infinite per cent on the ground. This irra- 
tionality is avoided by cutting off the 
spectrum at the eta energy E*, where 


pik, +2 07 orth Bev. (3.7) 
c Bets these assumptions, area ABCD in 
Es Eg a Fig. 2 has an effect of dy’ per cent on the 
Fig. 2 Variation of the differential intensity ground, where 
spectrum of primary cosmic rays by : } 
the electric field. E,: cut-off energy 1 = dy (%4)=——— I oe “dE. G3) 
at geomagnetic latitude 29, and E\ = ou 3 
By 1 Be. Then the decrease of area DOC! 
ee ee ere es | 
by the influence of sine 
Eg oleae where ats 
WA= 7 eae aE ae Ae je sa 
| Pa ares a tate a 


a Ea 69) extending ricinae 


‘ rval_ between Ey and By v 


| 
. 


=e 
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[B] In the case of I(E) dE=J-e-*" -dE 


(i) From Eq. (3.5), the decrease of cosmic ray intensity at the top of the 
atmosphere becomes 


Yo(%) = 100(1— e-*4*) + 100a *AAE-e*™) [—E,{—a(E,+4E)}] 


+100 a-4B-e "| — Se ee | (3.12) 
(E,+4E) | (E,44E)Y (E445 


where FE; denotes Logarithmic Integral, and 


A=2+2 (mpoc?)-a+5 (me)? FFL (mye?) B+ tetas 
2 3 
B=2 (mye?) +5 (mye? +11 a 23 (mye?) F Eb ee 


C=5 (mec’) - a (229C7)8- 


oe 
a sae i SM awatiahins 


D=11 (mc?) - Lees (myc*)*- 


a 
2\5 
3 +47 (myc?) B zi 3 +93 (myc?) Eigse 
(3. ‘ten 
(ii) The decrease of intensity on the ground:—By the same consideration as 


in the case of [A], the decrease of intensity becomes 


dy (%) = aes I (E+ 4E){1—L (E+ 4E)-4E} jw 


T(E) 
07 [, ME+4E)) ppp 0.7L (E+4E)-I(E+4E)-4E- 
ae [2 Te) | abe .. Te dE, (3.14) 


The first term on the right hand side of Eq. (3. dh is not pare to be=caleulatel 
under the above mentioned assumptions. For, ~ (= 


TW E+AE)=J-e° 4) = (J-e-94%)e-8#, : (3.15) 
Prem this équlition, we see that the spectrum of cosmic rays is not influenced by 


the presence of the electric field, and only its coefficient J changes to fe", Tine 
: other words, the intensity of the primary cosmic rays decreases without varying its 
_ intensity distribution at infinity, therefore this effect is eS at ay bo in 2 
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Y=t Ite (3.18) 
Figs. 3 and 4 represent the decrease of the vertical intensity in per cent, calculated 
from Eqs. (3.6), (3.10), (3.12) and (3.18) respectively. 

From these figures, we see that the latitude effect of the decrease of cosmic 
ray intensity is larger in the case of [A] than [B]. We may also conclude that the 
influence upon the decrease on the ground is more effective in the case of [B] than 
[Al 

In conclusion, we may expect the decrease of intensity of about one or two 
per cent on the ground by an electric field which makes cosmic rays lose their ener- 
gies of about 10° ev. On the other hand, in the case of the magnetic field, the varia- 
tion of the cut-off energy of about 10° ev. or more is necessary to produce such a 
decrease on the ground. So the decrease of cosmic ray intensity on the ground is 


more influenced by the electric field than by the magnetic field. 


3 
0 10 20 30 40 50 X 0 10 20 30 40 50 X 


Fig. 3 Decrease of vertical intensity on the Fig. 4 Decrease of vertical intensity at the 
ground in per cent, im~the case of upperzatmosphere, in per cent, in the 
AE=108 ev., as a function of geomag- case of AE=108 ey., as a function of 
netic latitude, for various energy geomagnetic latitude, for various ener- 
ie hehe gy spectrums. Curves I, II and III cor- 
ie Ae: ee eS mae respond to the Curves in Fig. 3, res- 

a=0.1 Bev.—, pectively. ; 


Curve III: I(E)dE=J-e—-*¥ .dE, 
, a=0.176 Bev.—, 
§ 4. The Magnetic Storm Model 


From the above calculation, the intensity of cosmic rays is decreased by the 
amount of about one or two per cent on the ground by the influence of the potential 
of about 10° volts. But we cannot adopt the electric field of the stream, adovocated 


a 
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by H. Alfvén, to explain the cosmic ray storm, because it is difficult to interpret 
experimental evidence (A) by this hypothesis. It is necessary to assume an electric 
field, growing with the current system. 

Of many theories of magnetic storm, it seems most appropriate to adopt that 
of S. Chapman and V.C.A. Ferraro.“ They suppose a neutral ionized stream, emitted 
from the sun with a velocity of about 2-10%cm/sec.. The earth entering into this 
stream, a westward current is produced. Its shape is uncertain, and for mathematical 
simplicity, it is said to be an equatorial ring current, spherical current or cylindrical 
current. 

It seems to be hopeless at present to determine more clearly the current system 
from the observations of the magnetic field alone. To examine the auroral phenomena 
during the magnetic storm@“.@%) is another way to determine the current system. 
Also, cosmic ray phenomena during the magnetic storm may be able to determine it. 

In the next, I will examine the current system from the point of view that the 
cosmic ray storm will be explained by the electric field. 

As pointed out by S. Chapman and V.C.A. Ferraro,“ or recently by D.F. 
Martyn,“® there must exist a radial electric field in the current system to stabilize 
itself. If an equatorial ring current is adopted, the potential difference, produced 

ws between the outer and the inner 

side of the ring current, is neces- 

sary to be more than 10° volts 
in order to produce the potential 
of about 10% volts on the earth. 
Then the energy of charged part- 


icle, incident upon the auroral 


zone from the ring??, is much 


| 
Bt 


larger than the observed one. But, 

if a closed current system, sur- 
rounding the earth, is adopted, 

x the above difficulty is avoided, 
- and furthermore the electric fied 


mts 


of this system influences most 
effectively. upon the cosmic rays. 
But as we have no knowledge 
about the current system, pro- 
duced practically, we must be 
satisfied with the next approxima- 
tion as follows: An infinite cy- 
lindrical current is adopted as an 


idealized current system for a 
Fig. 5 Cylindrical current’system as an idealized 


mathematical simplicity to esti- 
magnetic storm model. 


mate the electric field of ‘the 
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current, as shown in Fig. 5. We will then consider the effects produced by the dif- 
ference between the idealized system and the practical one. 

In conclusion, the current system, required from the stand-point of the electric 
field hypothesis of the cosmic rays storm, is that which produces a potential of the 


order of about 10% volts on the earth. 


§ 5. Estimation of Electric Field connected with Magnetic Storm 


Let us examine the electric field produced by the current system by using the 
above-mentioned idealized model. As pointed out in previous section, a radial electric 
field is required in the cylindrical current system for the stability of itself, as shown 
in Fig. 5. According to the storm theory of Chapman and Ferraro,“ the equations 
of the motion of charged particles in the current are 

—m,u,7/r = eK —eu,H, (5.1) 
—mue/r= —eH + eucH. (5.2) : 


m., m;: masses of electron and ion respectively, 
Ue, U;: velocities of electron and ion respectively, (west-ward positive), 
E: polarized electric field, 
H: magnetic field. 
Adding Eq. (5.1) and (6.2), 
: ; MuZ=M Ue +meue =reH (u;—Ue). (5.3) 
From Eq. (5.3), Meg > Mee , 
Moreover, rank Ge and Ferraro also deduce the next relation, 


(u U;—te) KU, OF Ue. , “i 


as 


‘The distance between the center ot the earth < the current system is E49 i 


i L ; 
> SIA Z=(8/Nmuz), ie — 64) 


SE Ogre Pe BEE ON | ee 
where — ; ee St SOBs > ; ® * i } - 
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If we neglect the magnetic field produced by the charged particles, and adopt appro- 
ximately the equation H=H,/Z*, we get the next relation from Eqs. (5.3); (5,4), (5.7) 
and (5.8). 


H.-4H-a  Z4 


dr — Hy, 
Ar (he ozs 
The electric field, produced in the current shell, is approximately 
5 ‘cP ) 
E=u,-H=u;,- oa (5.10) 
From Eqs. (5.9) and (5.10), the potential difference between A and B becomes (c.f. 
Fig. 5) 
V=E-dr=u;- PE eS 
4r B 


H,?-AH-a 
ae 


Ar 


oe 


_ (8/Nm,u,") : 


B ; 


Table I Calculated values of AE and Z 


Case (i) i) 


(Martyn’s estimation) 


(5.11) 


4; 
N AH(y) | us (cm/sec.) 1-108 2 - 108 3+ 108 4+ 108 

| AE(e.v.)= || 3.5-10° | 5:6-10¢ |.7.3-106 | 8.9- 108 

20 100 
ee 5.5 4.4 3.9 3.5 
yee or 0 vee cb aor thew, Ss 008 

ay a Te pena 

Coa 5.5 4.4 3.9 3.5 
AE(e.v.)= | 2.7+108 | 4.3-10¢ | 5.6+108 | 6.8- 106 

1 100 
iy Rge 4.2 3.4 2.9 2.7 


Case (ii) geeks. 


20 4x ° 


Ped 


(1/20 of Martyn’s estimation) 


AH(y) | us (cm/sec.) 1-108 2 - 108 
AE(e.v.) = 4.3 +107 
m ve Z= 3.4 
AE(e.v.) = 3.3 + 107 
2.6 


Case (iii) $= 50° 4x 


AH(y) 


L 


ux (cm/sec.) 


AE(e.v.) = 


AE(e.v.) = 
wh = 


Ai? 


1 - 108 


2-108 


(1/50 of Martyn’s estimation) 


3 - 108 


‘9.2.+ 107 
239 


FOxe LO% 
2.2 


1.5 + 108 
2.3 


il =)105 
1.8 


1.9 «108 


2.0 


1.5 - 108 
5 
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In this case, the potential, produced on the earth by this current system, equals to 
V in Eq. (5.11). Therefore, cosmic rays, arriving to the earth from infinity, lose their 


energies by an amount of 4E, respectively, where 
AE=eV. (5.12) 
Table I shows the calculated values of Eq. (5.12) for various values of 4H, N, u; 


and 8, where proton’s mass is adopted as the value of m;. ; 

But, it is for a mathematical simplicity and not essential to adopt such an 
infinite cylindrical current system. I will show this in the next. It is clear that the 
practical current system has no infinite length, but we have no knowledge about its 
shape or its dimensions. Therefore, by using an eqatorial ring current which is ex- 
tremely different in shape from the infinite cylindrical current, let us estimate the 
potential produced on the earth and compare with that of the cylindrical current. 
If the difference of both potentials will be small, it will be able to say that the 
potential produced on the earth by the practical current system will be not different 
very much from that of the infinite cylindrical current system. 

In the case of the ring current, the potential difference, produced between ‘the 


outer and the inner side of the ring current, is 


seopiabtige & “ts 
PRLS EES Pree MotEE™ (6.13) 
3 . 


In a case of N=20, 4H=10- Gauth, u,=2x 108 cm/sec., and B=H,?/4z, 
Vr=1.44 x 107 volts. (5.14) 
The potential produced on the earth by this ring current becomes approximately to 
one-sixth of Vr in Eq. (5.13) after some calculation. -Then 
: V=2.4 x 10° volts. - aa (5.15) 


From Table I, in the infinite cylindrical current, — 
V=5.6x 10% volts. — — @.16) 


Therefore, the potential produced on the earth in the case of the ring current equals — 
to about a half or one-third of that of the infinite cylindrical current. bay ie ae 
above calculation, it can-be said that the practical current system produces ft 
same ge! on the earth as the infinite Us Be current syste % 

planation of the Cosmic Ray Phenomena during the Ma 
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variations of mw; and p. 

It is undesirable to adopt such a large value of uw, as 2X 10° cm/sec., for its 
value is supposed to 2x10‘ cm/sec. from the magnetic storms. Decreasing of the 
vaiue of 8 is a way to increase the energy loss of cosmic ray, as shown in Table I. 
8 is a parameter for determining, to what distance from the earth’s center arrives 
the solar stream, and is given different values by Chapman and Martyn, respectively. 
If we can expect § to be one-twentieth or one-fiftieth of Martyn’s estimation, the 
decrease of the cosmic ray intensity will be explainable by such an electric field. But 
it may not be necessary for 8 to be ievested with such a small value because we are 
quite ignorant about the current system, so there may be something which makes 
the potential larger. 

(ii) As mentioned above, we cannot get a satisfactory result quantitatively, 
but it seems to be very reasonable to account for other phenomena by this electric 
field. 

(a) It seems to be suitab’e in explaining Exp. Evi. (C), because there are no 
cosmic rays arriving at any earth’s surface without suffering their energy losses by 
the electric field. This discussion holds true only if the potential of the earth is 
higher than that of infinity and does not depend upon the shape of the current sys- 
tem. But also in this case, the closed current system is most effective. 

(b) S.E. Forbush® pointed out that the decrease of cosmic ray intensity had 
latitude effect of several per cent on the ground. Theoretically, if J-E~* is adopted, 
the calculated latitude effect is too large, compared with the observed one. But if 
J-e-** is adopted, the calculated effect seems to be very reasonable to account for 
the observations, as in Fig. 3. : 

(c) As for Exp. Evi. (D), the interpretation is quite contrary to that of the 
magnetic field hypothesis,“ as stated in Sect. 1. From Eq. (5.11), the larger the 
radius of the current system, the smaller the decrease of cosmic ray intensity. But, 
we must not forget, the above expression only holds in the case of N=constant. 

(d) There is no explanation for Exp. Evi. (F) at present, but it does not seem 
' to me to be impossible to search for it. 

(e) Aithough the cause of the diurnal variation of cosmic rays is obscure at 
present, it may be sure that the diurnal variation is produced by the low energy 
region of cosmic rays. When the intensity of cosmic rays decreases with the growth 
of the electric field of the current system, the ratio of the low energy part to the 
total intensity decreases, i.e., the amplitude of the diurnal variation decreases slight- 
ly (cf. Fig. 2). Although this decrease will be far smaller than that of the variation 
of magnetic field, we cannot explain the increase of the amplitude of the diurnal 
variation of cosmic rays by the electric field of.the current system. It is necesssary 
to assume another cause to explain this phenomenon. 

; In conclusion, it seems to be very reasonable to account for the cosmic ray 
phenomena during the magnetic storm by such an electric field, except for the in- 
sufficiency of the decrease of intensity. ' é 
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§ 7. Observations which will determine the True Cause of Cosmic Storm 


As mentioned in a previous section, everything is working smoothly, in the case 
of the electric field hypothesis than the magnetic one. But we must not forget that 
the observations at the higher atmosphere during the magnetic storms play important 
roles, in judging the relative merits of these hypotheses. 

If the decrease of cosmic ray intensity should be explained by an improved 
Chapman’s magnetic field hypothesis, what should be the decrease at the upper 
atmosphere? ‘This is estimated from the latitude effect on the ground, as the varia- 
tion of magnetic field will cause only that of the cut-off energy. In Fig. 6, the de- 
creases of intensity at the upper atmos- 
phere are shown for various energy spect- 
rums and latitudes, when that of the 
ground is just one per cent. At the same 
time, Fig. 6 also contains that of the 
electric field hypothesis. But the latter is 
a rough approximation for we neglect the 
variation of the cut-off energy. 

From these figures, we see that the 
decrease of intensity at the upper atmos- 
phere is larger in the case of the magnetic 
field hypothesis than the electric one. If 
we are able to know the decreases at 
various altitudes from the observations, it 
will be decided which is the true cause. 

There is another way to decide the 
true cause. In the case of the magnetic 
field hypothesis, it is low energy region 
of cosmic rays that is prevented to arrive 
0 at the earth by the increase of the earth’s 


0 10 20 30 40 50 Xx magnetic dipole-moment, and therefore 


Fig. 7 Assumed decrease of vertical intensity there is no variation on the higher energy 
at the upper atmosphere, in per cent, 


: : ; region. 
as a function of geomagnetic latitude, ; 
for various energy spectrums, when On the other hand, in the case of 
that of the ground ‘s just one per the electric field hypothesis, the decrease 


cent. 

Curve I: J(E)dE=J-E-*-dE, 

a=2.75, of low energy but also in that of high 

Curve II: I(E)dE=J-e—«".dE, eae 2 
a=0,1 Bev! energy (c.f. Fig. 2). So the observation 

Curve Ill: J(£)dE=J-e—*".dE, : A ‘ 

«=0.176 Bev-1 Of high energy region of cosmic rays dur- 
and subscripts M and Erepresent the ing the magnetic storm is a key to 
curves, obtained from the magnetic 
field hypothesis and the electric one, ; 
respectively. Unfortunately, we cannot do it at 


of intensity occurs not only in a region 


determine the true cause. 
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once for the absence of such data, so it is most desirable to make such observations 
in future. 


§ 8. Conclusion 


It must be emphasized that adopting the infinite cylindrical current as the 
west-ward current system during the magnetic storm is only for a mathematical 
simplicity and is never essential. The discussions mentioned in the previous hold 
good also strictly in the case of the practical current system, because the electric 
field in the current system produces a potential on the earth, though less effective 
than the idealized one, as stated in Sect. 5. 

The difficulty with which this electric field hypothesis is confronted at present, 
is that the decrease of cosmic ray intensity is of the order ten or twenty times too 
large to be accounted for by the electric field, deduced from the Chapman’s storm 
theory. 

Except for this, it seems to be very reasonable to account for the cosmic ray 
phenomena, connected with the magnetic storm, by this electric field hypothesis. Thus, 
the relation between the magnetic and cosmic ray storm seems to be accounted for 
from the standpoint of the electric field, but the connection of the auroral phenomena 
with them remains obscure. 

Lastly, next two observations of the cosmic rays;—(i) the decrease of the 
high energy component of cosmic rays during the magnetic storm, and (ii) the de- 
creases of the cosmic rays at higher altitudes during the magnetic storm,—are most 
important for determining the true cause of the cosmic ray storm. It is most desira- 


ble to make such observations and determine the cause of the cosmic ray storm. 
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Letters to Editor 


On a self-exciting process in magneto-hydrodynamics 


Recenily, a simple and plausible explanation of the non-dipole part of the earth’s 
magnetic field and its secular variations was proposed by W. M. Elsasser and E. C. 
Bullard. The variations are attributed to hydrodynamical eddies in the earth’s core 
and to changes in electric current systems in the outer part of the core caused by 
the eddies. This success leads to the idea that the whole magnetic field of the earth 
is ultimately produced and maintained by a mechanism of induction. Elsasser dis- 
cusses the relation between the field and motion in a very general way and concludes 
that there may exist chains of relations in which a field H, interacts with a motion 
V, to produce a field H, which, in turn, interacts with another motion V, to reproduce 
H,. If it is possible to show the existence of such a “ self-exciting dynamo” in the 
earth’s core, the main field might also be explained. All that is necessary isa weak 
initial field to start with, and the field will then grow automatically until the energy 
dissipated as Joule’s heat becomes equal to that supplied. 

On the other hand, T. G. Cowling showed that the self-exciting process is 
impossible if the motion of the fluid is rotationally symmetrical about an axis and is 
confined to meridian planes. In addition, Elsasser showed that the same is also true 
for purely toroidal flow. Many believe that these theorems are foretastes of a stron- 
ger theorem that will prove the self-exciting process is impossible at all. In this 
connection, Bullard proposed a convection current model, which does not violate the 
above theorems, and by which the self-exciting dynamo may be possible. We have 
examined the problem whether the magneto-hydrodynamical equations have a solution 
representing the self-exciting process which is appropriate to the Bullard’s model. 

The problem is reduced to solve the following differential equations for the 
radial parts of the magnetic fields S,°, T,°, T,° and fie 

(16/3)-[d2S,°/dr?+-4/r-dS,°/dr|—(A8 x 24/5) -7° V2" To**=0 
(24/5)-[7-d°T!/dr? +6-d Ts! /dr|—(15/5)-Ane-dV;"/dr-S," 
— (288 x 6/35) +4 -7—"[ d(r! V2" T2”*)/dr + Ar V2"°)/dr +7 To?]=0, 
(48 x 6/5) -[7-d? T2*°/dr’? +6 -dT,*°/dr]+(96 x 6/5) 4x Vi°r Tr’* 
—(288 x 6/35) «4x -7—*[d(' Vo"° Tx) /dr + dr? V2?°)/dr-r T° =0, 
(48x 6/5) -[r-@2 Ts2*/dy? +6-dTy*/dr|—(96 x 6/5)-4ne- Vir Ti?” —(96 x 2/5)- 
Ank -1—[393 Vo" d(dS,°/dr+2/r-S\°)/dr + d{ dr? V2"°)/dr-Si°}/dr]=0, (1) 
where « is the electric conductivity of the earth’s core and 7 is the distance from the 
center of the earth. The radial parts of the fluid velocity fields V,° and V,* are put 


tentatively 
Vio (n= a Vio [1—(7/a)], V2(7) = a *V2.2d1—(7/a)P, (2) 


where a is the radius of the earth’s core and V;.o and V2.2. are undetermined cons- 
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tants of the dimensions of velocity. The boundary conditions imposed on the radial 
functions S,°, Tx®, T° and Tx’ are T:°(7)=T2*-(”)=0 (3), Si(7)+7/3-dS,°/dr=0 (4) 
at y=a. The radial functions should also be finite at 7=0. Thus our problem is 
reduced to solve (1) with (2) under the boundary conditions (3) and (4). It is 
easily seen that this is an eigen-value problem for Vi. and V2oe- If we get real 
eigen-values for Vi.» and V2.2. by solving the problem, we can thus prove that the 
process proposed by Bullard is self-exciting. An approximate method is used to obtain 
the eigen-values, and later this approximate method is shown to be good enough for 
studying the present problem. The eigen-values 4z*aV (where V!, denote the maxi- 
mum values of the velocities concerned) thus obtained are of the orders of 10° c.g.s. 
em.u. These results are to be compared with those obtained by Elsasser and 
Bullard in their discussions on the secular variation of the earth’s magnetic field, 
ie., <=3. 10 em.u., V~10—, a=3. 10°, 4nxaV~10*. The results obtained by us show 
that a self-exciting dynamo is possible by which the earth’s magnetic field is produced 
and maintained. The self-exciting dynamo is considered to have the similar structure 
as supposed by Elsasser and Bullared in their studies on the secular variation of the 
earth’s magnetic field. 


H. Takeuchi and Y. Shimazu 


Geophysical Institute, Faculty of Science, 
Tokyo University, Tokyo, Japan. 
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Drifting of Es Ionic Clouds 


Drifting of Es ionic clouds in Japan was already discussed by the author (1). 
It was obtained in the study that the direction of the drifting changed semi-diurnally 
with a mean speed of about 100m/s and that southward drifting during 19-0lh. 
was the most frequently and conspicuous. From an estimation for relations between 
the drifting and terrestrial magnetism, it was assumed in the report that the south- 
ward drifting during this period of 19-O0lh., in the northern hemisphere, might occur 
in region of higher latitudes than the center of an electric current that produced 
geomagnetic daily variation, and that northward drifting might occur in the lower 
latitudes during this period,—in the southern hemisphere, these relations would be 
reverse. 

In this report, some examples surporting this estimation are obtained from 


observational data in foreign stations besides Japanese stations. 


12-13 (FEB. 1950) 16 -17 


FEBRUARY 1950 


ise ee 341.20 0 6 12 18 0 UT 

Fig. 1 /Es variations representing northward Fig. 2 Variation-curves of monthly median 
drifting of Es at Brisbane, Canberra values of fEs representing northward 

and Hobart. drifting of Es during 19-Olh.. B, C and 

H of big lines are respectively Brisbane, 

Northward drifting of Es at three Canberra and Hobart situated on almost 

. oe eridi same meridian, and W and Ch of fine 
stations situated on almost same meridian, Fee cre Warhersginaird ACHR 
Brisbane (153.0°E, 27.5°S), Canberra situated in the neighbourhood of those 


(149.°E, 35.3°S) and Hobart (147.4°E, three stations. 
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42.8°S) in Australia, as shown in Fig. 1, may occur frequently during 19-01h. in summer. 

This tendency is seen in variation-curves of monthly median values of fEs too, 
at these stations and Watheroo (115.9°E, 30.3°S) and Christchurch (172.6°E, 435°); 
as shown in Fig. 2. Its mean speed in these examples is about 100-150m/s. As 
these stations seem to situate in higher latitudes than the center of the dynamo- 
current of the southern hemisphere in summer, such northward drifting is favourable 
to above-mentioned estimation. 

On the other hand, variation-curves of monthly median values of fEs at 
Wakkanai (141.7°E, 45.4°N), Fukaura (139.9°E, 40.6°N), Shibata (139.3°E, 37.9°N), Lan- 
chow (103.8°E, 36.1°N), Koku- 
bunji (139.5°E, 35.7°N), Yama- 
gawa (130.6°E, 31.2°N), Chung- 
Ta king (106.8°E, 29.4°N), Okinawa 
(127.8°E, 26.3°N), Guam (144.8°E, 
J ae 13.5°N) and Leyte (125.0°E, 

SA BT Brien K 11.0°N) are compared in Fig. 3, 
and their remarkable variations 

a during 19-Olh. extracted from 
the figure will be seen in Fig. 4. 
Southward drifting at first 


be eae seven stations is conspicuous in 


the figure as already reported, 
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and northward drifting seems 
to occur at Guam and Leyte 
with a speed of comparable 
Order abovementioned, though 
fine variations of drifting in 
each day at these two stations 
can not be obtained because of 
lack of values in each day. These 


Bo aCe 6 (5 ut 19 0 5 examples also may be favourable 
Fig.3 Variation-curves Fig. 4 Deviations from dotted t . * 
ee a: ot - : 
Be cuihly scotia lined ih baeh araceeshee Penks to the above-estimated relations 
values of fEs. Sta- of Fig. 3, concerning the In conclusion, the author 
tions are Wakkanai, geographic latitude of each : : . 
Get deep) Shibata, i. wishes to express his sincere 
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Errata 


K. MAEDA and M. IZAKI, Some Applications of Modern Mathematical 
Statistics to the Lunar Tidal Analysis of Ionospheric Height. 
(Vol. IN, No. 1, April, 1951) 
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